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APPLICATION OF REMOTELY SENSEI LAND-USE INFORMATION 
TO IMPROVE ESTIMATES OF STREW FLOW CHARACTERISTICS 


by Edward J. Pluhowski 
Abstract 

Land-use data derived from high-altitude photography and satellite 
imagery are presented for 49 basins in Delaware, and eastern Maryland and 
Virginia. Based on 1:100,000 scale maps from high-altitude photography, 
basin land cover was extracted at the generalized Level I and the more 
detailed Level II classification categories. Level I land-use data sum- 
maries were prepared for 46 of the basins using the 1:250,000 scale maps 
derived from Landsat imagery. Land cover in the basins ranged from 93.9 
percent urban at Little Falls Branch near Bethesda, Maryland, to 96.2 
percent agricultural at Morgan Creek near Kennedy vi lie, Maryland. 

Applying multiple regression techniques to a network of gaging stations 
monitoring runoff from 39 of the basins, it was demonstrated that land-use 
data from high-altitude photography provides an effective means of signif- 
icantly improving estimates of streamflow. Forty streamflow^characteri sties 
equations incorporating remotely sensed land-use information, were 
compared with a control set of equations using map derived land cover. 
Siginificant improvement was detected In six equations where Level I data 
was added and in five equations where Level II information was utilized. 

Only four equations were improved significantly using land-use data derived 
from Landsat imagery. Significant losses in accuracy due to the use of 
remotely sensed land-use information were detected only in estimates of 
flood peaks. Losses in accuracy for flood peaks were probably due to land 
cover changes associated with temporal differences among the primary land- 
use data sources. 
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INTRODUCTION 


Since 1888 when systematic streamflow records were first collected 
by the U.S. Geological Survey, more than 16,000 sites have been gaged 
in the United States (Carter and Davidian, 1968). Surface-water data are 
used for many purposes such as evaluating the water supply available to 
a town or city, designing bridges and culverts, or assessing the flood 
potential along a particular watercourse. A well designed stream-gaging 
network is of considerable value in studies attempting to assess man's 
impact on the hydrologic cycle. For example, urbanization will change 
streamflow patterns because of street paving, home and building con- 
struction, and the installation of storm sewers. These and other activities 
needed to develop urban environments alter important basin characteristics 
such as infiltration rates, generated, volume of storm flow, and the time 
required for water to move from any point in the basin to stream channels. 
Ideally, continuous streamflow monitoring v/ould be required before, during, 
and after development to appraise the impact of urbanization on, a .particu- 
lar watercourse. 

The general objective of the streamflow data program is to provide 
users v/ith water data at any site on any stream. Clearly, ft is neither 

i 

practical nor desirable to gage every site where data are required. It is, 
however, frequently possible to transfer streamflow information on un- 
regulated streams to other natural stream sites in areas of similar 
climatic and geologic settings. Thomas and Benson (1970) outlined a 
multiple-regression method of streamflow generalization. This procedure 
involves regressing a single streamflow characteristic (such as mean 
annual discharge) against the physiographic and climatologic characteristics 
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of gaged basins within a selected region. Equations obtained from the 
multiple-regression procedure contain only statistically significant 
basin characteristics, and the regression equations enable users to 
compute streamflow patterns at any site on natural streams within the 
region. 

Using basin characteristics derived from climatologic data and 
maps, detailed formulas were obtained by the multiple-regression 
procedure for a wide range of streamflow characteristics throughout the 
Nation, The results of these investigations, published in open-file 
reports, are available at the 46 district offices of the U.S. 

Geological Survey except Hawaii (Benson and Carter, 1973). The purpose 
of this investigation is to investigate the potential improvement of 
streamflow estimates by using land-use information obtained from 
high-altitude photographs and. satel lite images. Remotely sensed data 
to be tested were obtained from U.S. Geological Survey land-use maps 
compiled by the Central Atlantic Regional Ecological Test Site (CARETS) 
project. 

CARETS PROJECT 

The CARETS project was sponsored jointly by the National 
Aeronautics and Space Administration (NASA) and the U.S. Geological 
Survey. The principal objective of CARETS was to test the extent to 
which various remote sensor data systems could be used as input to a 
regional land-resources information data base (Alexander, 1974). The 
CARETS region covers 46,434 mi (74,712 km ) which includes Delaware, 
southern New Jersey, southeastern Pennsylvania, District of Columbia, 
and eastern Maryland and Virginia (fig. 1). 
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Figure 1, — Map of CARETS area showing location of basins for which 
land use was delineated. 




Figure 1. -- Map of CARETS area showing location of basins for which 
land use was delineated. 






NASA aircraft flown at altitudes of about 60,000 ft (18,300 m) 
provided colbr and color infrared photographs of the site in 1970 and 
again in 1972. The bulk of the high-altitude land-use analysis was done 
using the 1970 aerial photographs. However, parts of the site were 
masked by clouds in the 1970 high-altitude photographs and other aerial 
photographs taken as close as possible to the dates of the 1970 missions 
were required to complete land-use mapping of the site. Landsat-1 imagery 
was available at 18-day intervals following launching of the satellite 
in July 1972. Land-use mapping predicated on satellite imagery was derived 
from Landsat-1 data obtained principally during September and October 1972 
(K. Fitzpatrick, oral commun., 1976). 

Photo interpreters examined each piece of film or imagery for the 
major land-use types such as urban land, agricultural land, forests, wetlands, 
or water. Urban land is recognized by the patterns of buildings, houses, 
road networks, railroads, and other man-made features. The complex urban 
setting contrasts strongly on high-altitude photographs and images with the 
less complicated appearance of agricultural fields, forests, wetlands, and 
water. 

Land-use maps based on high-altitude photographs were produced at a 
scale of 1:100,000. Owing to resolution differences between Landsat 
imagery and high-altitude photographs, land-use maps derived from 
satellite imagery were prepared at a scale of 1:250,000. Forty eight 
sheets depicting land use of the CARETS area at a scale of 1:100,000 
and eight sheets at a scale of 1:250,000 have been released to- the 
U.S. Geological Survey open files, along with many additional map types 
to assist users in applying the data to land-use planning and environmental 
interpretation (Alexander and others, 1975). 
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LAND-USE CLASSIFICATION SYSTEM 


The classification system used in the CARETS project was one 
developed by a special interagency committee (R. Alexander, written 
commun., 1976) later slightly modified into the USGS Land-Use 
Classification System for use with remote-sensor data (Anderson and 
others, 1972). The scheme is a multilevel, hierarchical classification 
system which specifies the first two levels (table 1), and leaves the 
more detailed levels for later definition. Level I contains generaliied 
categories suitable for delineation from satellite imagery. Level II 
yields greater detail within each Level I category and is most suitably 
obtained using high-altitude photographs as a primary source. 


6 



Table 1 . — Land-use 

Level I Categories 
URBAN & BUILT-UP 


AGRICULTURAL 

FOREST LAND 
WATER 

NONFORESTED WETLAND 
BARREN LAND 


categories used in CARETS data base 


Level II Category Numbers and Titles 

11- Res1dent1al 

12- Commercial and services 

13- Industrial 

14- Extractive 

1 5- T ransportation , communi cati ons , 
and utilities 

16- Institutional 

17- Strip and clustered settlement 

18- Mlxed 

19- Open and other 

21 - Crop! and and pasture 

22- Orchards » groves, bush fruits, 
vineyards, and horticultural 
areas 

23- Feeding operations 

24- Other 

41 - Heavy crown cover (40% & over) 

42- Light crown cover (10% to 40%) 

51 - Streams and waterways 

52- Lakes 

53- Reservoirs 

54- Bays and estuaries 

55- Other 

61 - Vegetated 

62- Bare 

72- Sand other than beaches 

73- Bare exposed rock 

74- Beaches ^ 

75- Other 
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LAND USE IN SELECTED BASINS 


Using maps prepared in accordance with the CARETS classification 
system (table 1), land use was defined for selected basins listed in 
table 2. The basins for which land-use information is presented are in 
the northwest and north-central part of the CARETS region (fig. 1). They 
represent a broad spectrum of land cover ranging from predominantly 
agricultural in the Delmarva Peninsula to urban in the Washington-Baltimore- 
Wilmington corridor. Land-use data were obtained by drawing the boundaries 
of each selected basin on clear plastic sheets. These basin outlines, 
prepared at scales of 1:100,000 or 1:250,000, were used as overlays, on 
CARETS land-use maps. The percentage of a basin ascribed to any particular 
category was determined manually using a dot planimeter. The dot planimeter 
is a uniform grid of dots on a clear plastic sheet which was placed over 
the basin boundary overlay. Land use beneath each dot was recorded, the 
number of dots subtotaled by category, each category subtotal was then 
divided by the sum total of dots within the basin boundaries, and the 
result multiplied by 100 to yield percent. 

Land Use Based on High-Altitude Photographs 

Land-use information for 49 basins based on high-altitude photographs 
is summarized in tables 3 and 4 at Levels I and II respectively. At 
the 1:100,000 scale used to compile tables 3 and 4, the smallest 
depictable area is about 10 acres (4 hectares), or the equivalent of a 
square 656 ft (200 m) on a side (Alexander, 1975, written communication). 
Table 3, which shows generalized Level I land-use categories, is a 
compilation of the more detailed Level II category listings in table 4. 


8 



Table 2. -- Drainage basins analyzed for land use and gaging stations used in multiple regression analysis 


KO 



station 

Ko. 

Station name 

Latitude 

Longitude 

Drainage 

area 

(oi2) 

Period of record 
analyzed 
(water years) 

01477800 

it 

Shellpot Creek at Wilmington, Del. 

39®45'39" 

75°31’10‘' 

7.46 

1945-67 

01478000 

* 

Christina R. at Coocha Bridge Del. 

39°38’16" 

75®43'46" 

20.5 

1943-67 

01478500 

* 

Whi.te Clay Creek above Newark, Del. 

39°42'50" 

75°45'35" 

66.7 

1952-59, 1962-67 

01479000 

White Clay Creek nr. Newark, Del. 

39°42'00" 

75°41'10" 

87.8 

AA 

1952-56 ,1957-67 

01483200 

A 

Blackbird Creek at Blackbird, Del. 

39°21’58" 

75°40'10" 

3.85 

01483500 

A 

Leipsic River nr. Cheawold, Del. 

39°13'58" 

75°37’57" 

9.35 

1932-33,1943-57, 1958- 
67*’* 

01484300 

A 

Sowbridge Branch nr. Milton^ Del. 

38°48'51" 

75°19’39" 

7.08 

1956-67 

01484500 

A 

Stocklcy Branch nr. Stockley, Del. 

38°38'19" 

75®20'31" 

5.24 

1943-67 

01485000 

A 

Pocomoke River nr. Willards, Md. 

38°23'20" 

75°19'30" 

60.5 

1949-67 

01485500 

A 

Hassawango Creek nr. Snow Hill, Nd. 

38°13'45" 

75®28’20" 

44.9 

1949-67 

01486000 

A 

Kanokln Brook nr. Princess Ann, Md. 

38°12'50" 

75°40'18" 

5.8 

1951-67 

01486500 

V A 

Beavcrdam Creek nr. Salisbury, Md. 

38°21'05'' 

75°34'11" 

19.5 

1930-33,1934-35,1936-67 

01487000 

Nanticoke River nr. Bridgeville, Del.* 

38°43’42" 

75°33'44" 

75.4 

1943-67 

01487500 

A 

Trap Pond Outlet nr. Laurel, Del. 

38°31'40" 

75°29’00" 

16.7 

1951-67 

01488500 

A 

^!arshy hope Creek nr. Adamsvllle, Del. 

38°51'00" 

75°40'29" 

44.8 

1943-67 

01489000 

A 

Faulkner Br. at Federalsburg, Md. 

38®42'45" 

75°47’35" 

7.1 

1950-67 

01490000 

Chlcamacomico River nr* Saletn^ Hd. 

38® 30 '45" 

75®52'50" 

15.0 

1951-67 

01491000 

A 

Choptank River nr. Greensboro, Md. 

38®59'50" 

75®47’10" 

113 

1948-67 

01492000 

A 

Beaverdam Branch at Matthews, Md. 

38®48'40" 

75°58’15" 

5.85 

1950-67 

01492500 

Salle Harris Cr. nr. Carmichael, Md. 

38°57'55" 

76°06’30" 

8.09 

1951-56,1957-67** 

01493000 

A 

Unicom Branch nr. Millington, Md. 

39°15'00' 

75°51’40" 

22.3 

1948-67 

01493500 

A 

Morgan Creek nr. Kennedy ville , Md. 

39°16'50" 

76®00'55" 

10.5 

1951-67 

01494000 

Southeast Creek at Church Hill, Md. 

39°07'57" 

75®58'51" 

12.5 

1951-56,1957-65** 





Table 2. — Drainage basins analyzed for land use -- Continued 






Drainage 

Period of record 

Station 

Station name 

Latitude, 

Longitude 

area 

analyzed 

No. 




(tai^) 

(water years) 

01495000 

Big Elk Creek at Elk Mills, Md.* 

39°39'26" 

75°49'20" 

52.6 

1932-67 

01495500 

little Elk Creek at Childs, Md.* 

39°38’30" 

75°52’00" 

26.8 

1949-58 

01496000 

Nortneast Creek nr. Leslie, Md. 

39°37'40" 

75°56'40" 

24.3 


01579000 

Basin Run at Liberty Grove, Md. 

39°39'30" 

76°06'10" 

5.31 

1948-58,1965-67 

01586000 

N. Br. Patapsco R. at Cedarhurst, Md.* 

39°30’00" 

76°53'00" 

56.6' 

1945-67 

01589300 

Gwynos Falls at Villa Nova, Md. 

39°20’43" 

76'’44’01" 

32.5 

1957-67 

01590000 

North River nr. Annapolis, Md.* 

38°59’09" 

76°37'2-" 

8.5 

1932-67 

No gage 

Rhoae River nr. Galesville, Md. 

38°52'00" 

76°31'00" 

14.8 


01591000 

Patuxent River nr. Unity, Md.* " 

39°14'18" 

77°03'23" 

34.8 

1944-67 

01594500 

Western Branch nr. Largo, Md.* 

38°52'34" 

76°47'54" 

30.2 

1950-67 

01594600 

Cocktown Cr. nr. Huntington, Md.* 

38°38'27" 

76°38'07" 

3.85 

1957-67 

01594800 

St. Leonard Cr. nr. St. Leonard, Md. 

38°26’57" 

76°29'43" 

6.73 


01645200 

Watts Branch at Rockville, Md.* 

39°05‘03" 

77°10'38" 

3.70 

1957-67 

01646200 

Scott Run nr. McLean, Va. 

38°57'32" 

77°12'21" 

4.69 


01646550 

Little Falls Br. nr. Bethesda, Md.* 

38°57’27" 

77°06’31" 

4.1 

1944-59,1960-61, 1961-6; 

01648000 

Rock Cr. at Sherrill Dr., Washington, D.C.* 

38°58'21" 

77°02’25" 

62.2 

1928-67 

01649500 

N.E. Br. Anacostia R. at River dale, Md.* 

38"57'37" 

76°55'34" 

72.8 

1938-67 

01650500 

N.W. Br. Anacostia R. nr. Colesvllle, Md.* 

39°03’55" 

77°01'48" 

21,1 

1924-67 

01652610 

Holmes Run nr. Annandale, Va. 

38°50’47" 

77°10'2S’’ 

7.10 


01653500 

Henson Creek at Oxon Hill, Md.* 

38°47'05" 

76°58'42" 

• 16.7 

1948-67 

01653900 

Accotink Cr. nr, Fairfax, Va. 

38°48'46" 

77°13'43" 

23.5 


01655500 

Cedar Run nr. Warrenton, Va, 

38'’44'30" 

77°47'15" 

13.0 

1950-67 

01656800 

Cub Run nr. Chantilly, Va. 

38°S4'30" 

77°28'01’’ 

7.13 


01656940 

Cub Run at Lee Highway nr. Chantilly, Va. 

38^49 *59" 

77°27’50" 

39.6 










Table 2. -- Drainage basins analyzed for land use -- Continued 




Station Station name Latitude 

Ko. 

01657800 Giles Run nr. Woodbrldge, Va. 38°40’48" 

01658000 Mattawooan Cr. nr. Poiaonkey, Md. 38°35'45” 

A 

Station used in regression analyses. 

Annual maximum discharge only. 

^ Incluees entire drainage basin above confine ice with West 


Drainage Period of record 



analyzed 
(water years) 

1950-67 









Table 3. --Level I land-use classifications, in percent, for 
selected basins in Delaware, eastern .Maryland and 
Virginia, (Based on high-altitude photography). 


Index 

No. 

(fig. 2) 

STATION NAME 

LEVEL I Categories 

URBAN 


FOREST 



BARREN 

4778 

Shellpot Creek at Wilmington, 


1 






Del.-^ 

84.9 

3.5 

11.0 

0.6 

0 

D 

4780 

Christiana .River at Coochs 








Bridge, Cel.* 

20 . 9 

59 . 9 

19.2 

0 

0 

0 

4735 

White Clay Creek above Newark, 








Del. 

5.0 

78.0 

19.0 

0.05 

0 

D 

4790 

White Clay Creek nr. Newark,, Del. 

11.1 

69.7 

19.1 

0 

0 

0 

4832 

Blackbird Creek at Blackbird , 








Del.* 

0 

61 . 6 

37.6 

0.8 

0 


4835 

•k 

Leipsic River near Cheswold , Del, 

0 

82 . 3 

17 . 7 

0 

0 

0 

4843 

Sowbridge Branch near Milton, Del 

0 

46.9 

52 . 5 

0.6 

0 

0 

4845 

Stockley Branch at Stockley, Del.* 

1.3 

56.5 

42.2 

0 

0 

0 

4850 

k 

Pocomoke River near Willards, Md ■ 

0.2 

49.6 

50.2 

0 

0 

0 

4855 

Nassawango Creek near Snow Hill, 








Md.* 

0,2 

20 . 2 

79.6 

0 

0 

0 

4860 

Manokin Br. near Princess Ann,Md.“ 

0 

31 . 6 

68 . 4 

0 

0 

0 

4865 

Beaverdam Creek near Salisbury, 








Md.* 

5 . 1 

44.7 

49.8 

0.4 

0 

0 

4870 

Nanticoke River near Bridgeville, 








Del.* 

1 . 1 

57.6 

41 . 3 

0 

0 

0 


P O 













ORIGINAL PAGE IS 
DE POOR QUALITY 


Table 3. --Level I lajid-use classifications, in percent, for 
selected basins in Delaware, eastern Maryland and 
Virginia- -con tinned 


Index 

No. 

(fig. 2} 


LEVEL I Categories 

STATION NAME 

URBAN 


FOREST 

WATER - 

WETLAND 

BARREN 

4875 

Tran Pond Outlet near Laurel, 
Del. 

0 

26.3 

72.8 

0.6 

0.3 

Q 

4885 

Marshy Hope Ct near Adamsville, 
Del.* 

0 . 1 

58.0 

41.9 

0 

0 

0 

4890 

Faulkner Branch at Federalsburg , 
Md. 

0 

72,5 

27 . 5 

0 

0 

0 

4900 

Chicamaconico River, near 
Salem, Md,* 

0 

53.0 

46.8 

0 . 2 

0 

0 

4910 

Chdptank Ri^ er near Greensboro, 
Md.* 

0.1 

55.8 

43.9 

0 

0.2 

0 

4920 

Beaverdaia Branch at Matthews, 
Md.* 

0 

71.2 

28.8 

0 

0 

0 

4925 

Sallie Harris Cr.near Carmichael, 
Md.* 

0 

67.9 

32.1 

0 

0 

0 

4930 

Unicorn Branch near Millington, 
Md.* 

0.4 

70.1 

29.2 

0.3 

0 

0 

4935 

Morgan Creek near Kennedyvil le , 
Md.* 

0 

96.2 

3.8 

0 

0 

0 

4940 

Southeast Ct.at Church Hill, Md.* 

0 

73.5 

26 . 5 

0 

0 

0 

4950 

Big Elk' Creek at Elk Mills, Md.* 

1.1 

85.9 

13.0 

0 

0 

0 

4955 

Little Elk Cr. at Childs, Md.* 

1.2 

79.5 

18.3 

0 1 

0.9 

0 











Tsole 3. --levei I land-use classifications, in percent, for 
selected basins in Delaware, eastern Maryland and 
V irginia — Continued 


Index 
No . 

(fig. 2) 


LEVEL I Categories 

STATION NAME 

URBAN 


FOREST 

water 

WETLAND 

BARREN 

4960 

Nortneast Creek nr. Leslie, Md. 

4 . 0 

so . 1 

15.9 

0 

0 

0 

5790 

Basin Run at Liberty Gro\e, Md.* 

1 . 9 

73.4 

24.7 

0 

0 

0 

S 860 

North Branch Patapsco River at 
Cedarhurst, Md.* 

3.6 

70.9 

25.3 

0.1 

0 

0.1 

5893 

Gwynns Falls at Villa Nova, Md.* 

35 . 4 

23. 7 

40.0 

0 

0 

0.9 

5900 

North River near Annapolis, Md, 

0 

33.0 

67.0 

0 

0 

0 

5905 

* 

Rhode River nr. Galesville, Md. 

3 . 7 

39 . 7 

42.9 

12.2 

1.5 

0 

5910 

Patuxent River near Unity, Md,* 

1.5 

66 . 3 

32.2 

0 

0 

0 

5945 

Western Branch near Largo, Md. 

18.5 

38 . 6 

42.9 

0 

0 

0 

5946 

Cocktown Creek near Huntington, 
Md. 

1.6 

57.7 

40.7 

0 

0 

0 

S 946 

St , Leonard Creek near St . Leonard 
Md,* 

t 

0 

18 ' 9 

81 . 1 

0 

0 

0 

6452 

j5* 

Watts Branch at Rockville, Md. 

42,5 

40.9 

16.6 

0 

0 

0 

6462 

Scott Run near McLean, Va, 

55.6 

9 . 6 

34 . 8 

0 

0 

0 

64655 

Little Falls Branch near 
Bethesda, Md.* 

93.9 

0 

6.1 

0 

0 

0 

6480 

Rock Creek at Sherrill Drive, 
Washington, D.C.* 

53 . 3 

26.1 

20,3 

0,3 

0 

0 









Table 3. 


U1 


O O 
8 § 
I? 

Is 


--Level 1 land-use classifications, in percent, for 
selected basins in Delaware, eastern Maryland and 
Virginia-- Con tinned 


Index 
No . 



LEVEL I Categories 



(fig. 2) 

STATION NAME 

URBAN 


FOREST 

WATER 

WETLAND 

BARREN 

6495 

N.E. Br, Anacostia River at 
Riverdale, Md.* 

45.1 

IS . 8 

38 . 9 

0.2 

0 

0 

6505 

N . W Br. Anacostia River near 
Colesville, Md.* 

26.0 

42.4 

51.6 

0 

0 

0 

65261 

Holmes Run near Annandale, Va. 

67.9 

1 . 3 

30. 8 

0 

0 

0 

6535 

Henson Creek at Oxon Hill, Md.* 

63.2 

4.8 

52.0 

0 

0 

0 

6539 

Accotink Cr. near Fairfax, Va. 

71.1 

8.6 

20 . 2 

0 

0 

0 


Cedar Run near Warrenton, Va., 

1 .9 

63 . 1 

3 t . o 

0.4 

0 

0 

6568 

Cub Run near Chantilly, Va. 

49.9 

28.0 

22.1 

0 

0 

0 

65694 

Cub Run at Lee Highway near 
Chantilly, Va. 

18 . 5 

46.9 

34 . 0 

0 

0 

0 

6578 

Giles Run near Woodbridge, Va. 

11.8 

33.3 

54.9 

0 

0 

0 

6580 

Mattawoman Cr.nr. Pomonkey , Md. 

7.2 

24 7 

68.0 

0.1 

0 

0 


‘Station used in regression analy 
+Includes entire drainage basin a 

s 

ses 

bove con 

fluence v 

i th West 

River . 














Table 4. 


l,evel II land-use classifications, in percent, for selected basins in Delaware, eastern Maryland and Virginia 
(Based on high-altitude photography) . 


Category 

Number 

Level II 
Category 
Description 

c 
• o 

U fj 

u w«o 

O E 

Ci,r-« at 
iH »H 
1-4 ^ 

w Qi 
t/i ti Q 

Christiana R. 
at Coochs Br, 
Del. #4780 

White Clay Cr. 
above Newark i 
Del. #4785 

u 

u - 

X M o 
cO nS ct 
—< 

U 0) rj- 
'Z. 

I-* 

•H rt ^ 

2s c: Q 

Blackbird Cr. 
at Blackbird, 
Del, #4832 

i-4 

Kd 
• :s 

Oi W K> 

o 00 

W t3* 

•H ^ at 

u 

p. 

'H * ^ 

fl; J-» cj 
^ C Q 

u 

ca • 

C tn 
<l) O *a- 
OO ^ 00 

H *r4 3t 

f-4 S 

■§ . ■ 

w fid 

U X 
PQ O lA 

X-^ CO 
A) a ^ 

^ o ^ 

JA 4J 
t) to 

O r-4 

♦J a 

10 A O 

(0 
• -a 

os, u d 

A lO 

A) t-H DO 
^ r-C ^ 

O H te 

E 3= 
o 
o 
O 

cz 

1 

n 

Beaverdam Cr, 
nr. Salisbury, 
Hd. #4 865 

//// 

URBAN 

inm 

mmm 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn 

n 

Rcsidcntia 1 


IRSI!I 

1.8 

7. 1 






0.1 


1 1 

12 

Commercial 

2.9 

0.3 


0.8 









13 

liiilustri al 


2,0 

0.3 

1 . 5 



■* 






mm 

Extractive 


0,1 











IS 

Transportation 

1.5 

1 . 5 

0.1 

0.1 








3.2 

16 

Institutional 

5.4“ 

0.9 

0.2 

0.9 

~ 



1.3 





17 

Strip orclustered 


0.4 






1 

0.2 

0.1 


0.8 

18 

Mixed 

0.7 



0.3 









19 

Open or other 

4 6 

0.7 

0.3 

0.5 



“ 






mn 

AGRICULTURE 

II nn 

////// 

nnn 

nnn 

nnn 

nnii 

nnn 

WBBBm 

nun 

nnn 

nnn 

nnn 

21 


3.5 

59.9 


69.7 

61.6 

82.3 

46.5 

56 • S 

49.4 

19.2 

31.6 

44.3 

22 

Orchards 






■■■■ 

■■■■ 




• 


mm 



■■ 








1.0 


0.4 

Bl 






' 


0.4 


0 2 




//// 

FORESTLAND 

nun 

nun 

nnn 

nnn 

nnn 

nun 

nnn 

nnn 

nnn 

nnn 

nnn 

■nun 

41 


10. I 

18.2 

18.5 

18.7 

11 6 

17.5 

51,8 

MHWi 

49.3 

74.3 

67.9 

47.9 

ma 

Light crown cover 

0.9 

1.0 

0.5 

0.4 


0.2 

0.7 

!■■■ 

0 9 

5.3 

0.5 

1.9 

mm 

WATER 

KMHM 


moon 


wtum 

inura 

ilfJUl 


warn 

ihkhh 

innsn 

mm 

52 

Lakes 


mtnn^ 


mjmm 

mmm 








S3 

R.eservoirs 

0.6 


0.05 


0.8 

' 




- 


1 0.4 

54 

Bays and estuaries 













mm 

WETLANDS 

nun 

linn 

nnn 

nun 

linn 

nnn 

nun 

linn 

nnn 

nnn 

nnn 

\niin 

1 61 

Vegetated 




0 1 

’ 








iwm 

BARREN LAND 

nini 

nini 

mm 

mam 

mnm 

wnam 

wmm 

liMIM 

\mum 

wasm 

mmm 

WXKfllmtk\ 

1 75 

"Other" 




I^Hi 

!■■■ 

!■■■ 

flHHIli 

!■■■ 


\mmm 

\mam 

.HHH! 
































































































































































OF POOR QUAUTH 


Table 4. — Level II land-use classifications, in percent, for selected basins in Delaware, eastern Maryland and Virginia 

(Based on high-altitude photography). Continued, 


Category 

Nuffiocr 

Level II 
Category 
Description 

Nanticoke R. 
nr . Bridgevi 1 1 c 
Del. ff4870 

Trap Pond 
outlet nr. 
Laurel, 

Del. »4875 

u 

u 

O 

> 00 
o VI CO 
eC e 

X -o 

V) < 

^ a 

Faulkner Br, 
at Federal slitt-g 
Md. »4890 

o 

w a 

H C O 

e »-« o 
o cd cri 

U C/) ^ 

m % 

a • 

u e 

•rt * 

u e: s 

o 

• o 
o: X 
l/J o 
u: C ^ 
C U C) 
C V 
4^ k ac 

o • . 

DCS 

1 

n 

1 
M CO 
CQ C O 

C cn 

U ^ •ft 

O SEE 

o S! 

•»-i ♦ , 

fS -o 

3 c 2 ; 

Morgan Cr.nr. 
Kennedyvi lie 
Md. #4935 

Southeast Cr. 
nt Church Hill 
Md. M940 

Big Elk Cr. 
nt Elk Mills 
Md. M9S0 

UN 

URBAN 

nun 

unn 

nnn 

nnn 

////// 

JMWtMtm 


fMm§m 

WMMMtMM 

wum 

Tnrrr 

mm 

11 

Residential 

0.6 


■ 

' 


IHHi 



hhh 

HHH 


0.3 

12 

Coinnicrci a 1 

0.3 












mm 

Industrial 














Extractive 





HUH 

■HH 




hhh 


HHHi' 

HU 

_TjcajispQrtation 







||HH|| 

■hh 



HHOI 


LiAj 

Institutional 




■ 








0.2 

HQ 


0.4 


0. 1 



0. 1 



0.4 



0.6 

18 

Mixed 












||H^^ 

19 

Open or other 









■■■H 

|HH| 

HHH 


ma 

AGRICULTURE 

mmnm 

////// 

mam 

uiutm 

nnn 

nuKfl 

nnn 

nnn 

iniii 

nnn 

nnn 

mm 

IS 

Cropland 6 pasture 

mmm 

26. 1 


BPeWQ 

S3. 0 

55 . 3 

Hnm 

67 9 

70.1 

96. 2 

' 7 3 • S 

S5.9 

mm 

Orchards 

0.1 






HI^H 






23 














24 

Other 

0.1 

mBBi 












FORESTLAND 

nun 

nun 

nnn 

nnn 

mm 

wnnak 

nnn 

nun 

mm 


w^num 

WtlMMMtm 

m 

Heavy crown cover 

Hh|9 

IBBRI 

IHIKI 


44.0 

43.0 

28.8 

32.1 

27.7 



mmm 

' 42 

Light crown cover 

QiSEH 

9 7 

0.2 


2.8 

0.9 


' 

1.5 


0,2 

0.1 

imi 

IV’ATER 

wttnm 

nun 

nnn 

nnn 

mm 

nnn 

////// 

nnn 

mm 

mm 

mm 

mm 

52 

Lakes 





0.2 




0.3 




53 

Reservoirs 


0.6 












Bays and estuaries 













i!xm 

WETLANDS 

////// 

////// 

WUOiM 

mm 

mm 

nmi 

nnn 

nnn 

mm 

nnn 

nnn 

mm 

mm 

Vecetated 


0.3 




0. 2 


~ 





mm 

BARREN LAND 

mnm 

ntnm 

Wfaaank 

maaam 

nomoi 

waaaaai. 

waaaaatL 

wmtn 


mmm 

HUM 


Lll. 

"Other” 

‘ 







mm 

wnau^ 



rowwj 
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Table 4. — Level II Innd-uso rl,isslf1c.u Lons, in i>onotU, for st-ic. tod b.islns in Delaware, eastern Maryland and Vlrfiinxa 

(Based on high-altitude photography) — Continued. 


Category 

Number 

Level II 
Category 
Description 

WSti 

u 

u 

*J QJ O 
« »-< 

« ^ c» 
O w) 

jsi O Sk 
4-i i-S 

O Vi -O 
Z C JS 

Basin Run at 
Liberty Grove 
Md «5790 

O 

u in 

(/I v> 

& 3 

JZ G 

^ u o 

0 Id 00 
CU l/) 

• lU «c 
C-> 

. 

4-» *13 

:s cd s 

Ghynns Falls 
at Villa Nova 
Md. tfS893 

v« 

s 

v< 

> in o 
w i-i o 
Ci •“« cv 
o u* 
x: c: 

cc 

Vi C . 

o c *a 
z < 

o 
V. C-I 

4) fi 

> > Irt 
H V> O 
O Ct 

f-i m 
U OS 

*o o 
0 

x: ^ -c 
c 

CrJ 

A ° 
O 

c coo 
0) V< 

X rt at 
3 «-J 
<-> . 
fd Vi *3 
fr.es: 

1 

m 

— 

X o> 
u 

3 r-» 

h > 

00 ^ TT 
U vfr 
to O tt 
V* CC 
V* 

rt VJ ►xj 

3: « s: 

//// 

URBAN 

mill 

////// 

////// 

////// 

////// 

////// 

WiWJ 

igunm 

WMMMtKM 

'Mgm^ 

[nmg 

— 

11 


0.4 

0 5 

1.9 

1 7 

26.9 


E— 

1.0 

12,3 



26,2 

12 

Comraerci a 1 





2 8 


■■■■ 

0.1 

1 4 



7.2 

■ra 

Industrial 




B— 

0,1 









Extractive 




0.1 

0.3 




- 





Tr.ins aort at iOn 


0 7 



0.2 




0.5 



1 . 8 

riTi 

Institutional 





2 6 


0 3 


1,9 



5 0 

nra 

Strip or clustered 

0.4 

2 8 


0 6 

0.4 



0.4 

0 3 

1 6 



\mm 

Mixed 

0 4 



0 1 









I- 19 






2. 1 




2 1 



2.3 

mn 

AGRICULTURE 

imnni 

////// 

////// 

////// 

////// 

////// 

mill 

nnn 

nnn 

nnn 

nnn 

nnn 

L_21 

Cropland 6 pasture 

■^9.5 

80.1 

73.4 

70 9 

23 7 

33.0 

39.7 

66.3 

38,6 


IS. 9 

40.9 


Orchards 














Feeding operations 













BBBI 

Other 













mil 

FORESTLAND 

nun 

////// 

nun 

////// 

////// 

nun 

nnn 

nnn 

nnn 

WUMUk 

WMWtgMM 

rnttanm 

■a 

Heavy crown cover 

HBHl 

IE— 

22 7 

23.7 

33.3 

67.0 


31.8 

41.4 

40.7 

SO. 7 

11.8 

■a 

Light crown cover 

HRHB 

mam 

2.0 

1 6 

6 7 


'{■Ural 

'0.4 

1 5 


0.4 

4.8 

full 

WATER 

wnntk 

\nnn 

linn 

null 

linn 

////// 

linn 

nnn 

nnn 

nnn 

nnn 

nnn 

S2 



1^^ 

1^^ 


11^— 

1^— 




1^— 

IHBE 

IE— 

S3 

Reservoirs 

0.1 



0.1 




1 




1 

■ra 

Bays and estuaries 

bhhi 

IE— 

lEEH 

IE— 

IE— 

IE— 

1— 

[jniiiiiiiim 

i— — 

IE— 

11^— 

iiraiHn 

ran 

WETLANDS 

nun 

nun 

////// 

////// 

////// 

nun 

nnn 

nnn 

nnn 

nnn 

nnn 

nnn \ 

Rl 

Vegetated 

0.9 







IE— 

1— E 

IE— 

1— 

IE— 

ini 


!inm 


ii^ran 

lEHE 

IREE 

lEEl 


iR—ra 

inn— 

Iff— U 

If— 1 

If— 

Lz.s 

1 "Other" 

mmm 

jjraraE 

IE— 

0 1 

0 9 

^E— 

jQjjn 

IE— 

1— — 

IE— 

lEEi 

IBEE 
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Table 4. — Level II land-use classifications, in percent, for selected basins in Delaware, eastern Maryland and Virginia 

(Based on high-altitude photography). -- Continued. 


- >s 

h h 
0 4) 

«u g 
^ SJ 
C3Z 

o 

Level II 
Category 
Description 

Scott Run 
nr McLean 
Va. 06462 

Little Falls 
Br . nr 
Bethesda 
Md. #6465.5 

Rock Creek 
at Washington 
DC #6480 

C£, ^ 

<S 

*0 I/: 

U O'! 
k 4-» O 
02 (A > O 
O H 

♦ u e£ 

CU 

• c ^ ^ 

55 < 2 

4) 

• rH 

cs > tn 
• H /> O 
k 4-» G 1/V 

o o ^ 
o u 
ua ti 

' C M e-, 

V < c 3 


•— < 

• H 

k ::: tr# 
o to 
C 1/) 
C O vO 
o X u 
v> O 
c 

„ Cl 

Accotink Cr 
nr. Fairfax 
Va. #6 539 

c 

o 

C C lA 
r3 i> lA 
Oi k lA 

k o 
k «s 

G k « 

u c ^ 

>% 

w 

•w cO 
4-> nD 
C C tA 
3 « aO 

S4 ^ 

t) 

3k* 
U C ^ 

Cub Run 
nr Chantilly 
Va. #6569 4 

Giles Run 
‘nr Woodbridgc 
Va. #6578 

7/// 

URBAN 

mill 

mm 

mm 

mm 

mm 

Him 

wtnm 

faMtttm 

waMMtam 

mm 

wtMiaam 

ram 

mm 

Residential 

24 8 

HBHB 

■BEB9 

■m3 


■Em 

59.4 ' 



■■IIM 

■■■ 

H99I 

12 



12.0 

5.4 

HBm 

0 9 

mm 

5 4 

12.3 



0.3 



Industrial . 




HBSHI 


1.1 


1 8 





■n 

Extractive 



0. 1 

3 3 



6 2 



. 0 

0 7 

0.7 

wm 

Transportation 

12 6 


0 4 

2 


.3 . 5 

5 . 6 




7 . S 

1 . 4 

16 

Inst itutional 

0.7 


4.1 


0 6 

— 1 

6 9 

? ? 

2 0 

0.3 


0, 2 

.3 . 3 

17 

Strip or clustered 





0.4 


0.3 



0.6 

1 . 3 


13 

Mixed 



0,1 

mm 









19 

Open or other 

2.6 

9 0 

5 7 

mm 

9 2 

3.4 

4 . 1 

7.3 


34.9 

S. 2 


fim 

iCRTCUlTURE 

mum 

mm 

mm 

mm 

Him 

mm 


mm 

mm 

mm 

mm 

mm 


Cronland S pasture 

9 6 


25 7 

15 8 

42 4 

1 3 

4 8 

B 6 

63 1 

28.0 

46 9 

33 3 

Bis 

Orchard s 



IHRHI 




■■■ 

■■■ 


■■■■ 



23 

Feeding operations 













mu 
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Table 4. — Level II land-u<?e classlf icaclons, In percent, for selected basins in Delaware, eastern Jlaryland and Virginia 

(Based on hlgh^altitude photography), Continued, 
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For example, at Shellpot Creek (station No. 01477800), 84.9 percent of the 
basin is characterized by the Level I, URBAN category (table 3). This value 
was obtained by adding the various Level II categories listed under the 
generalized URBAN classification in table 4, Thus, the 84.9 percent URBAN 
Level I classification shown in table 3 for Shellpot Creek is equal to the 
sum of the following Level II URBAN categories listed in table 4: 


Category and number 

Percent 

Residential 

(11) 

69.8 

Commercial 

(12) 

2.9 

Transportation (15) 

1.5 

Institutional 

(16) 

5.4 

Mixed 

(18) 

0.7 

Open or other 

(19) 

4.6 


TOTAL 84.9 

Similarly, the Level I FOREST category for Shellpot Creek (11.0 percent) 
in table 3 was obtained by adding Level II forest heavy crown cover (10.1 
percent) and light crown cover (0.9 percent) in table 4. Because only single 
Level II categories, cropland and pasture, and reservoirs correspond to the 
general Level I category of AGRICULTURE and WATER respectively, identical 
values are shown at corresponding category levels for Shellpot Creek 
(tables 3 and 4). 

Based on high-altitude photographs, the highest measured percentage 
(93.9) Level I URBAN designation was at Little Falls Branch near 
Bethesda (table 3), By way of contrast, no urban development was detected 
in the high-altitude photographs of 11 Oelmarva Peninsula basins. Agricul- 
tural usage ranged from zero at Little Falls Branch to 96.2 percent at 
Morgan Creek near Kennedyville. Forest cover ranged from 3.8 percent at 
Morgan Creek to 81.1 percent at St. Leonard Creek near St. Leonard. With 
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the exception of Rhode River near Galesville, water areas identified in 
the 1:100,000 scale land-use maps amounted to less than 1 percent of the 
total drainage area of all basins. The Rhode River watershed is the only 
basin without a stream-gaging station as its downstream reference point. 
Land use given in tables 3 and 4 for the Rhode River catchment is for the 
entire basin above its confluence with West River. The high percentage 
(12.2) of the basin in the WATER category results from the largely es- 
tuarine lower part of the watershed. Wetlands were detected in four of 
the basins while only two basins had land use corresponding to the Level 
I BARREN category. 

Land Use Based on Landsat-1 Imagery 
The significantly lower resolution of Landsat imagery relative to 
high-altitude photography precludes its use as a data source for all 
Level II land-use categories. As previously noted, however, satellite 
imagery was used as the source base for preparing highly generalized 
Level I land-use maps at a scale of 1:250,000. The basic problem with 
Landsat imagery as used in this project is that Its spectral and tonal 
signatures cannot always be consistently matched with categories in 
land-use classification scheme^ especially where land parcels are small 
and categories are intermixed (Alexander, 1975, written communication). 
CARETS interpreters experienced particular difficulty in accurately 
mapping urban and built-up land in non-metropolitan areas using Landsat 
imagery (K. Fitzpatrick, 1976, oral communication). 
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Level I land use for 46 selected basins using satellite imagery as 
the primary source of land-cover information is shown in table 5. In 
general, these data are within 10 to 15 percent (by category) of the 
more accurate Level I land-use values based on high-altitude photography 
given in table 3. K. Fitzpatrick (written communication, 1975) reports 
that Level I land-use maps, when mapped from high-altitude photography 
were 7 percent more accurate for the entire CARETS area than the much 
less expensive Level I satellite-based land-use maps. However, accuracy 
differences greater than 7 percent between high-altitude and satellite 
sensors occur in table 5 owing partially to the small size of some of 
the basins selected for land-use analysis. Thus, in addition to errors 
stemming from lower resolution and problems with spectraT-signature 
discrimination, errors inherent in accurately positioning- such small 
basins on 1:250,000 scale land-use maps introduced additional variance, 
thereby further amplifying accuracy losses. Despite additional errors 
due to basin size, category differences in excess of 20 percent between 
Level I data based on high-altitude photography (table 3), and that based 
on satellite imagery (table 5) were detected in just eight basins. 

As anticipated, the largest discrepancies when comparing 
high-al titude with satellite sensor derived Level I categories generally 
occurred in suburban areas. Interpreters encountered difficulty segre- 
gating urban areas from surrounding non-urban land use in satellite 
imagery. For example, extensive urban areas in the N.W. Branch Anacostia 
River basin near Colesville, just north of Washington, D.C., were incor- 
rectly interpreted as agricultural land in Landsat-1 imagery;; accordingly, 
a high proportion of the basin (71 percent) was placed in the AGRICULTURE 
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Table 5. 


--Level I land-use classifications, in percent, for 
selected basins in Delahare, eastern Maryland and 
Virginia (Based on Landsat imagery) 


Index 

No. 

(fig-2) 

STATION NAME 

LEVEL I Categories 

URBAN 

AGRICUL- 

TURE 

FOREST 

WATER 

WETLA.NE 

4778 

Shellpot Creek at Wilmington, 
Del. 

86 

0 

14 

0 

0 

4780 

Christiana River at Coochs 
Bridge, Del.* 

12 

56 

30 


0 

4785 

White Clay Creek above Newark, 
Del.* 

0 

80 

20 

0 

« 1 

4832 

Blackbird Creek at Blackbird, 
Del.* 

0 

51 

49 

0 

[ 

0 

4835 

Leipsic River near Cheswold, Del, 

0 

91 

9 

0 

0 

4843 

Sowbndge Branch near Milton, Del.* 

0 

52 

48 

0 

0 

4845 

Stockley Branch at Stockley , Del. 

1 

61 

38 

0 

0 

4850 

Pocomoke River near Willards, 
Md,* 

0 

43 

56 

0 

1 

4855 

Nassawango Creek near Snow Hill, 
Md.* 

2 

26 

72 

0 

0 

4860 

Manokin Br . near Princess Ann, 
Md.* 

0 

30 

70 

0 

0 

4865 

Beaverdam Creek near Salisbury, 
Md.* 

4 

51 

45 

0 

0 

4870 

Nanticoke River near Bridgeville 
Del.* 

0 

54 

45 

1 

0 
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Table 5. 


-- Level I land-use classifications based on Landsat imagery, in percent — Continued. 


Index 



LEVEL I 

Categories 


1 

No. 

ffie. 21 

STATION NAME 

URBAN 

AGRICU- 

TURE 

FOREST 

WATER 

WETLAND 1 

4825 

Trap Pond Outlet near Laurel, Oel. 

0 

28 

72 

0 

1 

0 

4885 

A 

Marshy Hope Cr. near Adainsville , Del* 

4 

56 

40 

0 

0 

4890 

A 

Faulkner Branch at Federalsburg, Md. 

0 

71 

29 

0 

0 

4900 

A 

Chicamacomico River near Salem, Md, 

0 

47 

53 

0 

0 

4910 

A 

Choptank River near Greensboro, Md. 

1 

55 

44 

0 

0 

4920 

A 

Beaverdam Branch at Matthews, Md. 

0 

93 

7 

0 

0 

4925 

A 

Sallie Harris Cr. near Carmichael, Md. 

0 

69 

31 

0 

0 

4930 

A 

Unicom Branch near Millington, Md. 

0 

74" 

26 

0 

0 

4935 

A 

Morgan Creek near Kennedyville, Md. 

0 

97 

• 3 

0 

0 

4940 

A 

Southeast Cr. at Church Hill, Md* 

0 

71 

29 

0 

0 

4950 

A 

Big Elk Creek at Elk Mills, Md. 

0 

80 

20 

0 

0 

4955 

Little Elk Cr. at Childs, Md.* 

0 

62 

33 

0 

0 

4960 

A 

Northeast Creek nr. Leslie, Md. 

0 

71 

29 

0 

0 

5790 

A 

Basin Run at Liberty Grove, Md. 

0 

100 

0 

0 

0 i 

5860 

North Branch Patapsco River at 
Cedarhurst, Md. 

2 

81 

17 

0 

1 

“ 1 

5893 

A 

Gwynns Falls at Villa Nova, Md. 

45 

25 

30 

0 

0 I 
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Table 5. Level I land -use classifications based on Landsat inagery, in percent — Continued 


Index 

No. 

(£ig. 2) 
5900 
5910 

5945 

5946 
I 5948 

j 

6452 

6462 

64655 

6480 



LEVEL I Categories 


65261 

6535 

6539 

6555 

6568 


Kortn River near Annapolis, Md. 

* 

Patuxent River near Unity, 

•k 

5i?estern Branch near Largo, Md. 

k 

Cocktown Creek near Huntington, Md. 

St. Leonard Creek near St. Leonard, Md. 

k 

iCatts Branch at Rockville, Md. 

Scott Run near McLean, Va. 

k 

Little Falls Branch near Bethesda, Md. 

Rock Creek at Sherrill Drive, 
Washington, D.C. 


E. Br. Anacostia River at Riverdale.Mdl 55 


p.W. Br. Anacostia River near Colesville 
Md.* 

Holmes Run near Annandale, Va. 

Henson Creek at Oxon Hill, Md. 

Accotink Cr. near Fairfax, Va. 

Cedar Run near Warrenton, Va. 

Cub Run near Chantilly, Va. 












Table 5. — Level I lana-use classifications based on Landsat imagery, in percent — Continued 



Index 



LEVEL 1 

Cateeories 




No, 

ff2.2. 2^ 

STATION NAME 

URBAN 

AGRICUL- 

TURE 

FOREST 

WATER 

WETLAND 


6578 

Giles Run near Woodbridge, Va. 

49 

0 

51 

0 

0 


6580 

Kattawoman Cr. near Pomonkey, Md. 

1 

29 

70 

0 

0 



*Station used in regression analyses 
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category (table 5). Based on high-altitude photographs only 42 percent 
of the basin was agricultural and 26 percent was designated urban 
(table 3). Using land-use maps derived from satellite imagery only 
6 percent of the basin was categorized as urban (table 5). 
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EXPERIMENT DESIGN 


The approach used in evaluating remotely sensed land-use data as a 
means of improving streamflow estimates was based on (!) selecting as many 
stream-gaging stations from the basins listed in tables 3-5 as possible 
to perform a meaningful multiple-regression analysis, (2) applying the 
same basin and climatic characteristics utilized in the streamflow program 
analysis of the Maryland district of the U.S. Geological Survey (Forrest 
and Walker, 1970) to the study basins in order to develop regional equations 
needed to compute specific streamflow characteristics, (3) incorporating 
selected' Level I and Level II land-use categories developed from both high- 
altitude and satellite sensors to define other sets of streamflow equations, 
and (4) comparing standard errors of estimate for each streamflow characteristic 
(control) equation developed using the basin characteristics available to 
the Maryland district of the U.S: Geological Survey with those generated 
by incorporation of remotely sensed land- use information. 

STUDY BASINS 

Records of 10 or more years are generally required to develop 
meaningful streamflow statistics. Streamflow records spanning at least 
10 years were available for 39 of the 49 basins for which land-use infor- 
mation is presented (tables 3-5). These stations (table 2) formed the 
network of study basins selected for multiple-regression analysis. The 
study basins drain into the Chesapeake Bay, Delaware Bay and the Atlantic 
Ocean (fig. 2), and are situated in the Piedmont and Coastal Plain 
physiographic provinces. The boundary between these provinces trends 
northeast through the Washington-Baltimore-Wilmington urban corridor. The 
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Figure 2. — Map showing drainage patterns and the location of gaging stations analyzed 
in this report. 



Figure 2. — Map showing drainage patterns and the location of gaging stations analyzed 
in this report. 



Piedmont Is characterized by rolling topography, low hills and ridges, 
and fairly steep side slopes. The Coastal Plain is low, flat, and poorly 
drained on the Delmarva Peninsula, but west of the Chesapeake Bay is more 
rolling with slightly improved drainage. 

Average annual basinwide precipitation Is quite uniform throughout 
the area with the lowest amount of 39.9 inches (1010 mm) reported at 
Cedar Run near Warrenton, Va. and the highest amount of 47.0 inches 
(1190 mm) at three Delmarva Peninsula basins (table A-1 , col. 19). As 
previously noted, the study basins exhibit a wide variety of land cover 
ranging from primarily urban in the Washington, Baltimore, and Wilmington 
metropolitan areas, to extensively forested west of the Chesapeake Bay 
in the abandoned farm areas just beyond the limits of urban development, 
and to agricultural in much of the Delmarva Peninsula. 

STREAMFLOW CHARACTERISTICS 

The streamflow characteristics (dependent variable) used in the 
streamflow analysis of the Maryland district span the full range of 
discharge regimen observed at 105 gaging stations. These include 
measures of high and low flows, discharge variability, and long-term 
average monthly and annual streamflow. Forty streamflow characteristics, 
in cubic feet per second, evaluated using all or some of the 39 gaging 
stations In this report, are as follows; 

f 

Q mean annual discharge, defined as the arithmetic 

average of the annual mean flows. 

q mean monthly discharge, where the subscript refers 

to the numerical order of the month beginning with 
January as 1 , 
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standard deviations of the annual means, 





standard deviations of the monthly means, where the 
subscript n refers to the numerical order of the 
month beginning with January as 1, 
annual flood peak discharge at T-year recurrence 
interval; recurrence intervals of 2, 5, 10, 25, 
and 50 years are denoted as P5 P-|0 ^25 
PgQ respectively. 


flood volume characteristics are the annual highest 
average flow for 3-day periods at recurrence intervals 
of 2 and 25 years (V^ 2 25^’ 7-day periods 

at recurrence intervals of 2, 10, and 25 years (V^ « 



10 , ^ 7 , 25 ^’ 


Mq -j- low-flow characteristics are the annual minimum 7-day 
average flows at recurrence intervals of 2, 10, and 20 

years ”7.10, ”7,20>- 

Dgg discharge equaled or exceeded 50 percent of the time. 


BASIN CHARACTERISTICS 

Characteristics Based on Maps and Weather Records 


Correlation studies performed on the Maryland district streamflow 
analysis incorporated 12 independent physiographic and climatic parameters 
into the multiple regression analysis as follows: 

A, drainage area, in square miles, as shown in the latest 
U.S. Geological Survey streamflow reports, 

S, main-channel slope, in feet per mile, computed by the 
10- to 85-percent method (Benson, 1962), 

L, main-channel length, in miles, measured from gaging 
station to basin divide, 

E, mean basin elevation, in feet above mean sea level, 
measured from topographic maps by the grid method 
(Benson, 1962). 

OEIGIIJAL 
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L, main-channel length, in miles, measured from gaging 
station to basin divide, 

E, mean basin elevation, in feet above mean sea level, 
measured from topographic maps by the grid method 
(Martins, 1968), 

area of lakes, ponds, and swamps, in percent of total 
drainage area, determined by planimetering such areas 
on topographic, maps, 

F, forest area, in percent of total drainage area, 
measured from topographic m'ps by the grid method, 
soil index, a measure of po*.entia 1 maximum infiltration 
capacity, in inches, estima'.ed from data provided by the 
U.S. Soil Conservation Service, 

P, mean annual precipitation, in inches, determined from 
isoheytal maps prepared from National Weather Service 
records , 

^24,2, precipitation intensity, expected once every two years 
over 24-hour periods, in inches, estimated from U.S. 

Weather Bureau Technical Paper 29, 

S^ mean annual snowfall, in inches, from snowfall maps 
prepared from National Weather Service records, 
average minimum January temperature, in degrees Fahrenheit, 
from National Weather Service records, 

T 7 average minimum July temperature, in degrees Fahrenheit, 
from National Weather Service records. 


Characteristics Based on High-Altitude Photograph 


Land-use classifications based- on high-altitude aerial photograph 
were tested as independent variables in the multiple regression analysis. 
These classifications, expressed in percent of total drainage area, 
are as follows: 
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U,, Level I urban or built-up land which comprise areas 
of intensive use with much of the land covered by 
structures , 

U. Levels I and II agricultural land consisting 

Q y 

predominantly of croplands and pasture. 

Level I forested land, 

T 9 

U, Level I and II water areas includes total area covered 
w, 

by lakes, reservoirs, streams, and estuaries, 

U Level II, residential, consisting of housing ranging 

I » 

from high density (multiple-family units) to low 
density (houses on large lots), 

U, Level II, industrial, consisting of land devoted to 

* y 

light to heavy manufacturing. 

Level II, other urban or built-up land consisting of 
parks, cemetaries, zoos, waste dumps, golf courses, 
and undeveloped land within, an urban setting. 

U^-j Level II, forest land, light crown cover (10 to 40 
percent) , and 

Level II, forest land, heavy crown cover (40 percent 
or greater). 


Characteristics Based on Landsat imagery 


Land-use classifications based on Landsat-1 imagery were also tested 
as independent variables in the multiple regression analysis. These 
classifications, expressed in percent of total drainage area, are as 
follows: 


Z Level I urban or built-up land which comprise areas of 

W y 

intensive use with much of the land covered by structures, 
1 Level I agricultural land consisting predominantly of 
croplands and pasture, and 
Z- Level I forested land. 
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REGRESSION ANALYSIS 


The multiple regression technique used defines the relation between 
a single streamflow characteristic (dependent variable) and an array of 
climatic, physiographic, and land-use characteristics (independent variables) 
for a selected network of stream-gaging stations. Only those independent 
variables that account for significant measures of variance in the stream- 
flow characteristic under analysis are included in the regression equation. 
Those independent variables that had at least a 95-percent probability of 
effectiveness were deemed -significant to the equation. An indication of 
accuracy provided by the equation relating a streamfloW' characteristic to 
significant basin characteristics is provided by the standard error of 
estimate. The standard error of estimate is a range of error such that the 
value estimated by the regression equation is within this range at about two 
out of three sites, and is within twice this range at about 19 out of 20 
sites for the sample population. 

Stepforward multiple regression analyses were performed by digital 
computer using STATPAC program D0094. The program eliminated doubtful 
dependent variable entries, added a small constant (0.0001) to those 
dependent variables which go to zero, and transformed all dependent 
variables and selected independent variables to their logarithms. The 
independent variable that accounts for most of the variance in the depen- 
dent variable was identified and entered into the regression equation. 

Then the next most effective variable was added to the equation. Because 
the significance of an independent variable in the equation changes with 
the addition of each new variable, all previously included variables were 
retested with the addition of a new variable, and any variable shown to be 
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no longer significant was deleted from the equation. The addition of 
variables accounting for a progressively smaller part of the variance 
in the dependent variable continues until the equation is not signifi- 
cantly improved by the inclusion of any additional variables. For each 
streamflow characteristic equation, the program provided the multiple 

I 

correlation coefficient, percent of total sums of squares of the dependent 
variable that are explained by the regression, and the standard error of 
estimate of the dependent variable. Program D0094 also tabulated observed, 
computed, and residual values of all streamflow characteristics at each of 
the 39 gaging stations used in the analysis. 

Observed, calculated, and measured values of all dependent and inde- 
pendent variables used in the Maryland district streamflow analysis were 
obtained from the Streamflow/Basin Characteristics retrieval program E796 
and are listed for each station in table A-1 (cols. 1-7, 19-55, 57-66). 


MULTIPLE REGRESSION MODEL 

The model equation used in the multiple regression analyses is: 
log Y = b-| log + b2 log Xg + b^ log 

* ® * '>n+l^ n+1 * *’n+2* n+2 •••• 

or its equivalent form: 

Y = ^1 Xg *^2 X^ ^ Vl Vl 


^ ^n+2 ^n+2 


+ b X ] , 
m m 


where 

Y 


b, to 
1 m 


a streamflow characteristic 
basin characteristics 
regression constant, and 
regression coefficients. 
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In this analysis, X-| through were logarithmically transformed 
whereas through were not transformed prior to calculations. 
Independent variables which tend to vary widely, such as (A) drainage 
area and (L) main channel length, were log (base 10) transformed, where- 
as those subject to relatively small variations were used directly. In 
addition to drainage area and main channel length, (S) main channel 
slope and (E) mean basin elevation were log transformed. All other 
basin characteristics were relatively stable and were used directly in 
the model equation. The model equation was applied uniformly for the 
development of control and experimental equations without comparing its 
effectiveness as a predictive tool with models wherein all variables 
are logarithmically transformed. Rather, simple comparative tests were 
performed to evaluate the usefulness of remotely sensed land-use data 
in improving estimates of individual streamflow characteristics. 

Specifically, the model was applied to 39 gaged basins in the 
CARETS region where land-use maps based on high-altitude photography 
and satellite imagery are available. A control set of equations was 
developed using the same basin characteristics that Forrest and Walker 
(1970) incorporated into their evaluation of the Maryland district 
streamflow program. The regression model was then applied success- 
ively to each of three experiments where additional land-use data 
were incorporated as follows: 


OF POOR 
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(1) four level I land-use categories derived from 
high-altitude photography, 

(2) six individual and combined Level II land-use 
categories derived from high-altitude photography, 
and , 

(3) three Level I land-use categories derived from 
Landsat-I imagery. 

Comparisons were then made between the control equations for 
individual streamflow characteristics and those developed for each 
of the above experiments to determine whether significant improvement 
in the standard error of estimate had resulted in any of the 40 
streamflow-characteristic equations. Changes of 10 or more percent 
in the standard errors of estimates between the control and experi- 
mental equations were arbitrarily deemed significant. 

The remotely sensed land-use categories selected for analysis 
depended on frequency of occurrence and percent basinwide coverage 
of each category, and category accuracy relative to map derived 
land-use data. For example, only four of six possible Level I land-use 
categories based on high-altitude photography were tested in experiment 
1. The Level I categories of wetlands and barren land were not used 
because of the 39 basins in the regression analysis, wetlands were 
detected in only three basins and barren lands in just two basins 
(table 3). Moreover, with the exception of the Rhode River basin which 
was not used in the regression analysis, the portion of either category 
(wetlands or barren land) relative to total area in any of the basins 
was less than one percent (table 3). Map derived percentages of areas 
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covered by lakes, ponds, and streams were used in experiment 3 rather 
than remotely sensed water data based on satellite imagery. 

Resolution problems as v^ell as spectral and tonal signature degradati 
precluded accurate detection of the small water bodies found in most 
of the test basins. 

REGRESSION EQUATIONS 

Tables 6-9 summarize the results of the multiple regression 
analyses. The first column of each table indicates streamflow 
characteristic (Y) coded in accordance with the scheme developed on 
p. 31- The last column lists the regression constant (a) corres- 
ponding to a particular streamflow characteristic. Regression 
coefficients (b^) for those independent variables found to be 
significant at the 95-percent level are listed in the intervening 
columns. Not all 39 stations in the test network were used in 
defining each of the regression equations shown in tables 6-9. 

Owing to varying periods of operation and special purpose gages, 
sufficient data to define streamflow frequency relationships for all 
40 characteristics was not available at all gaging stations. For 
example, two of the gages were designed to measure floods 
(crest-gage stations) and were used only in the flood-peak compu- 
tations. The number of stations used to develop each streamflow 
characteristic equation is as follows: 
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ORIGINAL 
OF POOR QUALITY 



Streamflow Characteristic No« of stations 


^ 2 , ^5, ^0, ^25 ■ 

39 

^A, ^1-12, ^^1-12 

37 

•^7,2, ^^7,10 

34 

^7,10 

33 

*^7,20 

32 

°50 

29 

^3 >2 

26 

'^7,25 

25 

^3,25 

24 

^50 

15 


The regression analysis results incorporating physiographic and 
climatic basin characteristics identical to those used in the Maryland 
district analysis are listed in table 6. These are the control 
equations with which equations using remotely sensed land-use infor- 
mation were compared. Tables 7 and 8 present equations based on the 
inclusion of four Level I and six Level II land-use categories, res- 
pectively. These categories were based on land-use maps using 
high-altitude photographs as the primary Information source. Level I 
land-use data based on Landsat-1 imagery at three category levels were 
also analyzed and the results are listed in table 9. 
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Table 6. — Control equations obtained by regressing streamflow characteristics against physiographic and climatic basin parameters 

obtained £rom climatologlc data and USCS topographic maps. 
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Cr«>naEe Hsln 
iroa chana 


I A 


1.006 
1.027 
1.029 
1.022 
1.006 
,959 
.986 
. aS2 
1.022 
. 887 
.987 
,994 
1.038 


Main Main Mean 

channel channel basin 

slope length elevation 

5 I E 



forest; 

cover 

F 

Sci 2 
ir.ae»c 

Sa 

Mean 

annual 

precip 

r - 

Precip. 

intensity 

Irt<al 

Snowfal 1 


0,120 


0.108 



. 140 

0.019 




,096 

.017- 



.0017 

. 159 





. 112 


.161 





. 173 

0 .013 





.027 


' 



.021 


. 168 





JanJarv July Regross- 

ain.Kua maximum >or 

tenp temp conseane 

T, Ty a 









Table 6. Control equations obtained _bj^ , regressing streamflow characteristics against physiographic and cliaiatlc basin parameters 
obtained fron cllmatologlc data and USGS topographic maps*' — Continued. 
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Table 6. - obtained by regressing streamflow characteristics against pnysiographac and cllr^tlc basin parameters 

obtaxned from climatologlc data and l/SGS topographic maps Continued. 
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Eiperiaeiital equations obtained by regressing streamClow characteristics against physiographic ana ciiiaatio oasin parai»et®rB 
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Table 7- — ^xperljnental equations obtained by regressing streamlfow characteristics against physlographxc and climatic basin paraneters 

and four level I land-use categories derived from climatologlc data, USGS topographic maps, and high-altitude photographs 
— Continued. 
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Table 7. 




Model 


Experiinectal equations obtained by regressing streamflow characteristics agaltst physiographic and climatic basin parameters, 
and four level I land use-categories derived from climatologic data, USGS topographic maps, and high-altiUude photographs 
— Continued, 
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*Althou8h (January mean minimua temperature) was used in the regression analysis it was not significant; accordingly it was not 
Hated in this table 



Table 8#- — Experutental equations obtained by regressing streanflow characteristics against physiographic and clioatic paraneters, 

and Six level II land^use categories derived fron climatologic data, USGS topographic aaps, and high-^altitude photographs 




T*JJie 8. Experinental equations obtained by regressing streamflow characteristics against physiographic and climatic parareters, 

and six level II land-use categories derived from climatologic data, USGS topograpnc naps, ard hig'-altit^de protographs— Contm..e3 
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Tables , — Eaperuaental equations obtained by regressing streooClou characteristics against physiographic and climatic para-teters, 

* . and SIX level II land-use categories derived froa clinatologic data, dSGS topographic caps, a, id hign-altitaoc pfccrograpns -- Continued. 
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'Although Tj (January mean minimum temperature) WHS. ^Oaed in the regression analysis it was not signiticint, Uvcordingly, it was not listed 
in this table. 
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Table 9 , — Experlsental equations obtained by regressing streamflow characteristics against physiographic and cliiaatic basin paraiwters, 

,and three level I land-use categories 'derived froTi climatologic data, USGS topographic maps, and Landsat imagery. 
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Table 9. — Experimental equations obtained by regressing streamflow characteristics against physiographic and climatic basin para-reters, 
and three level I land-^se categories derived from climatologic data, USGS topographic maps, and Landsat imagery — coicimcd. 

W „i>z , Is bs b b 7 P-f^ bs 1. 4- Sn -r hfO~TJ ^ ba /7 -t bf2. i O/^ Di-i- 

Kodel IS Y ~ /-i S L E /O 


Exponent or coefficient of basin characteristic* 


Flow 

Drainage 

area 

Main 

channel 

slope 

Ham 

channel 

length 

Soil 

index 

Mean 

annual 

precip, 

Precip, 

intensity 

Snowfall 

July 

roaximiia 

temperature 

1 Level I Categories 


charact- , 
eristics j 

Urban 

Agricul- 
ture j 

Forest 

Regression 

constant 

y 

A 

S 

L 


P 

I (2V,z) 

Sn 

~7 

^u 

2a 

if 

£ 

SD, 1 

1 .022 



.217 


- 






-1.226 

SO, 1 

1 074 

1 

, i 

.194 








- 901 

SD^ 

1 .087 






- 023 





,200 

SO3 

1.037 




.032 


-.015 





-1 278 

SC^ 

1.018 




.021 ‘ 



j 




-1.076 

S05 

1,015 











- 223 

S05 

1.039 






i 

1 

1 





- .400 

SD, 

.922 










- DCSC 

- 040 

SOg 

1.080 











- 0S3 

SDg 

.858 


1 

i 




I 

.092 i 




-8,054 

'■’lo 

1,076 



-.281 



-.029 

! 

1 




.994 


1.031 


1 









- 313 

SD12 

1 1.087 


1 

! 







j 


- ,316 


equations obtained by regressing streaoiflow characteristics against pnysiographic and ctimatic basin paraaeters, 
and three level I land-use categories derived from climatologic data, USG3 topographic maps, and Undsat imagery — concjuiucd. 
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To illustrate the use of the regression equations, assume that the 
2-year peak flow (P^) is required for Shellpot Creek at Wilmington 
(fig. 2, index No. 4778} using (1) map and climate data (table 6), (2) 
added Level II land use based on high-altitude photography (table 8), 
or (3) added Level I land use from satellite imagery (table 9). The 
equations for (1) are: 
from table 6: 

= /\1*067 5O.77O (1.312 - 0.0089F - 0.023S^) 

from table A1 : 

p2 = 7 46^ 67 10 ■ 0-0089(19) - 0.023 (20)] 

p2 = 8.535 (25.50) 10°'®^^ 

P2 = 1050 ft^/s, 

( 2 ) 

from table 8: 

p ^ ^0.699 ^0.453 (1.067 + 0.0042 - 0.0065 

from table A1 : 

P2 = 7 64°'^^^ 271°-'^^^ 10 0.0042 (69.8) - 0.0065 (10.1 )] 

P2 = 4.143 (12.65) 10^-^^^ 

P2 = 1030 ft^/s, 
and (3) 
from table 9: 

p ^ ^0.991 3O.745 .jQ (2.022 - 0.397 I24 2) 
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from tabl e A1 : 

= 7 . 46 °-^^'' 67 . 1 °-^^® 10 
P^ = 7.326 (23.05) 

P^ = 870 ft^/s 

Each of these 2-year peak flow estimates at Shell pot Creek, based 
on regression analyses, is below the 1,200 ft^/s computed from actual 
station records (table A1 , col. 24). Part of the variation between 
predicted and recorded discharge is due to chance. However, Shellpot 
Creek drains a highly urban area and is subject to flash flooding owing 
to the impervious nature of its basin. Because the regression analysis 
is based on rural as well as urban streams, fairly sizeable discrepan- 
cies in the 2-year recurrence flood between actual and estimated values 
were anticipated at the station. 
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ACCURACY COMPARISONS 


Tables 10-12 identify the significant independent variables in both 
the control and experimental equation arrays as well as the standard error 
of each equation in logarithmic units and approximate equivalent percentage. 

The percentages represent arithmetic averages of the plus and minus percent 
of the mean, calculated using the standard error in log units. Thus, an 
average standard error of 18.5 percent, corresponding to 0.08 log units, 
represents a range of 20.2 percent on’ the plus (high) side and 16.8 on the 
minus (low) side of the streamflow characteristic mean (Hardison, 1969). The 
last two columns show the percent change in the standard error resulting from 
inclusion of land-use information in the analysis. Changes of 10 or more 
percent in the standard error of estimate are considered to be significant. 

Plus percent changes are indicative of improved accuracy whereas minus changes 
represent a loss of accuracy. Percent change values are given for all stream- 
flow characteristics except the three 7-day low-fl,ow categories. Less than 
50 percent of the variance in each of these categories was explained' by any 
of the 7-day low-flow regression equations. This strongly suggests that other 
unidentified independent variables should have been included In the regression 
analyses. Accordingly, conclusions regarding relative accuracy improvements were 
not made for any of the low-flow categories. 

Experiment 1 

In the first experimental array of regression equations, four of six 
possible Level I land-use classifications derived from high-altitude photo- 
graphs were tested; namely, FORESTLAND (U^), AGRICULTURAL (U^), URBAN AND BUILTUP 
(U^), and WATER (U^). As previously noted, the remaining two Level I categories. 
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Table 10 . — Comparison of standard error of estltoate changes resulting from inclusion in the regression analysis of 
four level I land-use categories derived from hlgh-altltude photography. 
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— CoDparlson of standard error of estimate charges resulting froia Incluslo’-. In the regression analysis of 
four level I land-use categories derived from hlgh-altltude pnotography — Continued, 
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Control eQuatlons | Experimental equations 


Standard error of estimate 


A, S„ 

A. S„, P 
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A, Pfc, p, Si 




• Ty . 
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in log units 
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Table 10. — Comparison of standard error of estimate changes resulting ^ 

four level I land-use categories derived from hlgh-altltude photography 


Flos 1 
haract- 
ristlcs 


Significant predictive variables 


Control equations | Experimen tal equations 


Standard error of estimate 


in log units 


Controx 


•In percent 


Control 


?ercer. 
change 

?lus I Hirus 


F; S; 8^ 
F; S, 


F; B, S^, T^ 


A, Ei Uf} Uy 
A; S, S^, ^2l4,S‘ 
A, Si S„. Uj.; U, 
A, Si 
A; P 


A} S^j Uf 


1,509 1.509 


No meaningful equations derived 
Vo meaningful equations derived 
No meaningful equations derived 


^ZU,2’ 
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Table u 


— Comparison of standard error estimate changes resulting from Inclusion In the regression analysis of 
six level II land -use categories derived from hlgh-altitude photography. 


Flovf 

charact- 

eristics 

Y 

Significant predictive variables 

Standard error of estimate 


Perce"t 

in log units 

i 

1 

in oercent 

cheir 


Control equations 

Experimental equations 

Control 

Exper . 
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Exoer. 
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A, Si; l2J),2 

A. Sl5 -2ll.2i 'Jfn 

0.062 

0.057 


lH.li 

13.2 


8.3 


Qi 

A; Si; ? 

A; Si, P , 

.061 
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i 

IH.I 

1.1.1. 


o 

0 

°-2 

A; Si; P 

A; Si; ? 

1 .067 

.067 


15.5 

}5.5 


Q 

& 


A; Si, F 

Ai Si; Ufh 

.073 

.072 

1 1 

16.9 

16.7 


1.2 


q!) 

A; Si; 

A; Si; l2H^2 

.072 

.072 


16.7 

16.7 


0 

0 


Ai Sjj; l2(|^2 

A; Sn; ^2^,2 

.095 

.095 


22.0 

22.0 


0 

0 

^6 

A 

A; E; Si; Ufh 

.133 

.100 


31.1 
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i 

25.4 


qj 

A i S„ 

A; Si; '^fh 
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33.5 • 
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<38 

A 

A 
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.120 
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0 

0 
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^n 
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1 

1 
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'*10 

A 

A 
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.139 
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0 

9 

‘ill 

A 

, A 
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.117 


27.3 • 
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0 

0 

^12 1 

A; S. 

i 

Ai Si ^ 
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.081 
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18.7 

I 

M.. 1 

0 
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TaMe 11 . — Comparison of atandard error estimate changes resulting from inclusion in the regression analysis of 
six level II land-use categories derived from hlgh-altltude photography — Continued. 
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A 
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■ A 
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Percent 
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Table 11. — Comparison of standard error estimate changes resulting from inclusion in the regression analysis of 
six level II land-use categories derived from blgh-<-altitude photography — Continued. 
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Percent 

change 

in log units 

1 
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Control equations 

Experimental equations 
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^2 
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^5 

A} F; S; 

n 

A; s. s„. 
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38 8 



10.2 

5*10 

A, F; 

A. s; 5^, 

,ll)7 

,17 A 
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41.1 



19.1 

^25 

A; F; S 

A. S; 

.158 j 
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37.3 

44.4 


1 

x7.i 

^50 

A, P 
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1 
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63.2 



0 

^3,2 . 

A 

A; Up 

i 

.126 

.106 
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24.6 


16.3 


''^3,25; 

A 

' A; Up 

.lil6 1 

.135 


3‘*.^ 

31.6 

1 
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-^7,2 

A 

A 
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.069 
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20.7 


0 

0 

^7,10 

A 

* A 

.103 

.103 
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23.9 


; 0 

0 

^7.25 

A 

A 

.102 
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23.7 


C 

i ^ 

”7,2 

A 

- A; Sj^; Ufjj 

.T91 

.685 


No meanlr 

1 

igful equation derlve< 

J 
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L 

L 

1.39*) 
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Ko neaningful equat 

Ion derived 
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A 

A 
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Table 12 Comparison of standard error of estimate changes resulting from inclusion in the regression analysis of tiuee level I la'-d-use 

categories derived from Landoat imagery. — Continued. 
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charact- 
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in log units 

1 

in percent 


:’£e 

Y 

Control equations 

Experimental equations 

Control 

Exper. 


Control 

Exocr . 


Plus 


SD 

a 

A. 

A; 

0.085 

0.085 


19.7 

19 7 


0 

0 

"°i 

A. Sj^ 

A. S, 

.09I1 

.09i4 


21,8 

21.8 


0 

0 


A; 

A.a„ 

.107 

.107 

■ 

2*1.9 

2*1 9 


0 

0 

SD^ 

‘ A. P 

1 

A. S„, P : 

.095 

.095 


22.0 

22 0 i 


0 

0 

SD^ 

A, P 

A. P 

,080 

.080 


16.5 

18. 5 


0 - 

0 

“5 

A 

A 

.116 

.116 


27.0 

27 0 


0 

0 

SDg 

A 

A 

.1I19 

.IU9 


35.0 

35.0 


0 i 

0 

SD^ 

A, F 

A, 

.19? 

.190 



i 1(5 1 


2.8 


SDg 

A 

A 

.153 

..153 


35.9 

. 35 9 


0 

0 

SDg 

A; • 

A, T^ 

.155 

•155 


36. li 

36.I1 


0 

0 

S“io 

A; S^, F. 

A. S,. 

.15a 

.157 


35.7 

36 9 



3 ** 

“11 

A 

A 

.15!j 

.15lt 


36.2 

36.2 


0 

0 

®^12 

A , 

A 

.137 

.137 


32.1 

32 1 


0 
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Companson of standard error of estimate changes resulting from inclusion in the regression analysis of tnree level I lana-use 
categories derived from Landsat imagery -- Continued 


Flow 

charact- 

eristics 

Y 

Significant 

predictive variables 

Standard error of estimate i 


Percent 
c 1 1 'I »” ■* 

in log units ! 

1 

] 

in percent j 

Control equations 

Experimental equations 

Control 

Exper. 

Control 

1 

Zw'oer , 

Pius 

r* 

L 


A, F. S, 

^24,2 

0.158 

0.188 


37.1 

44.6 

i 


20.2 


A, F, S, Sjj 

A; S, , 

.150 

.190 


35 2 

45.1 

1 

i 


22 1 

*”10 

A 5 S, 

A, 3, 

.147 

.198 


34.5 

47.1 



36 5 

‘’25 

A, r, s 1 

Aj 3, 

.158 

.191 


37 8 

45 4 

' 


20 1 


A.P 1 

A, P 

.259 

.259 


63 2 

63.2 


0 

0 

'^3,2 

A 

Zu 

.126 

.109 


29.4 

25.4 


13.6 


^3.25 

A 

A5Z„ 

.146 

.128 


34.2 

29 9 


12 6 


"t.2 I 

. A 

• 

.089 

.08h 

i 

■ 

20.7 

19.5 


5.3 


^.10 ' 

A 

A 

.103 

.103 


23 9 

23.9 


0 

0 

^T,25 

A 

A, Z 

.102 

.087 1 


23.7 

! 20 2 


l4,8 


*'V,2 

A 

A, 3 , Z 
* n’ u 

.791 

.651 


Ho raea-.l’igful e.iuation oci’ived 



^.iO 

t 

L 

L 
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1.394 


Ho meaningful equation deriveo 



" 7,25 

A 

A. 

1.509 

1.435 


Ho meonlncrtil c*iaation derivea 



"50 

A; 

^24,2’ 

' .143 

.119 

! 

33.5 

27 8 


17.0 




BARREN LAND and WETLAND, were identif.ed in only five of the 39 basins 
used in the correlation network, and were not included in the regression 
analysis. Throughout the analyses, was substituted for the USGS topo- 
graphic map (scale 1:24,000) derived forest (F) category which was used 
in the control equations, and was used in place of the USGS map derived 
storage (S^) category also used in the control equations. No substitutions 
were required for U or U because neither category was available for use 

Cl U , 

in the original (control) equations. 

Results of the experiment are listed in table 10, which shows that 

t 

11 equations were improved (six significantly) and five equations sus- 
tained a loss of accuracy (two significantly). By far the most often used 
independent variable in the regression analysis was FORESTLAND (U^) as 
indicated below: 

Streamflow Number Number of times that indicated 


Characteristic 

of 


variable occurred 


Type 

equations 

^a 

Uf 

% 

^w 

High 

10 

0 

4 

3 

2 

Average 

14 

1 

6 

0 

2 

Low 

3 

0 

1 

0 

0 

Variability 

13 

1 

1 

0 

1 

All characteristics 

40 

2 

12 

3 

5 


Five of the six streamflow characteristic equations significantly 
improved by inclusion of Level I land-use information involved mean flow 
characteristics (pg q^ q.|^ and Dgg) whereas one flood volume charac- 

teristic (V^ g) equation was similarly improved. A significant accuracy 
loss was detected in two flood peak characteristics (Pg ' Examin- 
ation of the four significant variables affecting the Pg relationships 
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(table 10) in both the control and experimental equation arrays indicates 
the presence of three dissimilar variables. By way of contrast, the three 
significant variables governing the P25 flood characteristic were identi- 
fied in both tests; however, F was used in the control set whereas was 
used in the experimental set. Owing to the loss of accuracy due to the 
inclusion of in the analysis, F (map derived) is the preferred inde- 
pendent variable for estimating 25-year flood peaks rather than which 
was obtained from high-altitude aircraft photography. 

Experiment 2 

In this experiment six Level II categories were included in the 
regression analysis to evaluate the possible impact of more detailed 
land-use information on streamflow estimates. As in Experiment 1, Level II 
data were derived from high-altitude photographs of the CARETS region. 

Two forest categories were included to depict heavy crown cover (U^^) and 
light crown cover (U^-j). Categories denoting residential (U^), industrial 
(Uj) and, open and other (U^) urban development were also incorporated in 
the analyses. The urban open and other (U^) category consists of golf 
courses, some parks, cemeteries, and undeveloped land within an urban 

I 

setting (Anderson and others, 1972). The last Level II classification 
used in the analysis was a combined cropland and pasture category (U^\ 
which essentially corresponded to the Level I agriculture category used 
in Experiment 1. Level II was not substituted for (map derived 
storage) in Experiment 2 because it appears that the category, based 
on 1:24,000 scale maps, portrays surface-water area with an equivalent 
accuracy to that derived from high-altitude photographs. 
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Thirteen equations were improved (five significantly) and five were 
reduced in accuracy (four significantly) by the inclusion of Level II 
land-use data derived from high-altitude photography (table 11). The 
independent variable most often appearing in the test equations was “fh 
whereas Uj never proved to be significant in any of the 40 equations as 


shown below: 
Streamflow 

Number 

Number of times that 

i ndi cated 

Characteristic 

type 

of 

equations 

U 

a 

variable occurred 
“fh “fl “r “o 


High 

10 

0 

4 0 3 

0 

0 

Average 

14 

1 

5 0 0 

0 

0 

Low 

3 

0 

1 0 0 

0 

0 

Variability 

13 

1 

1 3 1 

1 

0 

All characteristics 

40 

2 

113 4 

1 

0 

Not surprisingly. 

the results 

( 

of this 

test closely parallel 

those 


in Experiment 1 in that the streamflow characteristic equations signifi- 
cantly improved in Experiment 2 were identical to five of the six 
characteristics similarly improved in Experiment 1. Significant accuracy 
losses were sustained in four of the five flood peak characteristic 
equations as evidenced by large minus percent changes (10 to 19 percent) 
in the standard errors for these experimental equations. As in Experiment 
1, more accurate flood estimates were generated in the control equations 
where F (map derived forest cover) appears as a stronger independent 
variable than either or (aircraft derived forest categories). The 
use of Level II aircraft derived land use generated a slight overall loss 
in accuracy in the equations when compared with the Level I categories 
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used in Experiment 1. Thus, the use of more detailed land-use discrim- 
ination provided by Level II was unwarranted in this particular stream- 
gaging network. 

The loss of accuracy in estimating flood peak discharges at all 
frequency intervals except the 50-year return period, where forest cover 
is relatively unimportant, is probably a function of how well the land- 
use information represents the selected streamflow study period. For 
example, flood flow records used in this analysis included all available 
gaging-station records through September 30, 1967. The maps available 
for determining forest cover (F) in the control equations were prepared 
predominantly during the late 1950's which approximates the median period 
of actual data collection at the gaging stations (table 2). Land-use maps 
which were derived from high-altitude photographs obtained in 1970 and 
1972 reflect conditions beyond the streamflow analysis cutoff date. 

Because flood flows are highly dependent on forest cover, the values for 
this factor (F) used in the control equations were better suited as flood 
flow predictors than either the Level I (U^) or Level II and aircr 
aircraft derived forest cover estimates obtained three to five years 
beyond the flood analysis cutoff date. 

Experiment 3 

Owing to a significant loss of land-use detail in Landsat imagery, 
only three of six possible Level I categories were tested in Experiment 3. 
These include agriculture (Z^), forestland (Z^), and urban and built-up 
(Z ). As in Experiment 2, S (map derived storage) was retained to reflect 
the percentage of each basin covered by lakes, ponds, and swamps. Level I 
forestland (Z^) was substituted for map derived forest (Fj. and Z^ 
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represent land-use characteristics which were not considered in the control 
equations. Aside from the substitution of and the addition of and 
to the analysis, all other basin characteristics tested in Experiment 3 
were identical with those used in the control equations. 

Only six equations were improved (four significantly) and an identical 
number were reduced in accuracy (table 12). The independent variable 
appearing most often in the analyses was which was significant in a 
total of six low- and high-water equations as indicated below; 


Streamf 1 ow 

Number 

Number 'of times 

that indicated 

Characteristic 

of 


variable 

occurred 

type 

equations 

^a 

' ^f 

^u 

High 

10 

0 

2 

4 

Average 

14 

1 

0 

0 

Low 

3 

0 

0 

2 

Vari abi 1 i ty 

13 

0 

1 

0 

All characteristics 

40 

1 

o 

o 

6 


As in the high-altitude photography experiments, flood peak equations 
were adversely affected by inclusion of remotely sensed land-use infor- 
mation. Four of the five flood peak equations showed significant accuracy 
losses. Satellite forest cover, (Z^) obtained principally in late 1972, 
was not as effective as map derived values (F) in portraying conditions 
representative of the flood flow data analyzed in this report. Moreover, 
additional difficulties in land-use discrimination in satellite imagery 
that were not encountered in high-altitude photography introduced further 
errors in evaluating Z^. The combination of these and other error factors 
interacted to amplify flood-flow accuracy losses to a range of 20 to 36.5 
percent (table 12). 

original PAGE K 

OF POOR 



Using a network of gaged basins in the Delmarva Peninsula, Hollyday 
(1976) found that 12 streamflow characteristics were significantly 
improved with the inclusion of Landsat derived land-use information. 

Hollyday extracted the following categories from satellite imagery for 
use in a multiple regression analysis (1) forest, (2) riparian (streambank) 
vegetation, (3) water, and (4) combined agricultural and urban land use. 

Only one accuracy loss (December mean discharge) was detected in his 
regression analysis of 20 gaging stations, all of which were included in 
this study. 

SUMMARY AND CONCLUSIONS 

Maps incorporating the CARETS land-use classification system were 
utilized to determine land cover in selected basins of Delaware, eastern 
Maryland and Virginia. Land-use maps based on high-altitude photographs 
were used to prepare Level I (generalized) and Level II (more detailed) 
classifications for 49 basins. Only Level I classifications could be 
defined on the 1:250,000 scale maps derived from Landsat-1 images. Land 
use varied from highly urbanized in many basins in the Washington-Baltimore- 
Wilmington corridor to heavily agricultural in the Delmarva Peninsula. 

Using a network of gaging stations consisting of 39 of the 49 basins 
for which land cover was defined, it was demonstrated that land-use data 
derived from high-altitude aircraft photograohs are effective in signifi- 
cantly improving streamflow estimates. Significant improvement in 
accuracy, defined as a 10 or greater percentage reduction in the 
standard error of estimate, was detected by comparing streamflow 
characteristic "control" equations with three experimental equation sets. 

The control equation set consisted of basin characteristics used in a 
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review of the streamflow program of t.ie Maryland district of the 
U.S. Geological Survey. Land-use data based on high-altitude photographs 
and satellite imagery were used in the experimental equation sets. 
Comparisons of the experimental and control equations utilizing land-use 
information derived from high-altitude photographs showed significant 
improvement in six equations incorporating Level I data and in five 
equations where Level II categories were used. Only four equations 
showed significant improvement using land-use information derived from 
Landsat-1 imagery. The lower resolution of imagery relative to high- 
altitude photographs and difficulties in classifying certain spectral 
signatures tend to lower the effectiveness of satellite sensors as a 
means of providing detailed land-use information. 

Of the wide range of streamflow characteristics tested, remotely 
sensed land-use data yielded losses in accuracy only in estimates of 
flood peaks. These losses in accuracy were probably due to land cover 
changes stemming from temporal differences among the three primary land- 
use data sources. For example, high-altitude photographs and satellite 
imagery were obtained primarily in 1970 and 1972, respectively, and 
streamflow records analyzed in this study terminated on September 30, 

1967. Thus, remotely-sensed land-use data were not synchronous with 
the period of flood-flow analysis. By way of contrast, map derived 
land-use data incorporated in the control equations were obtained primarily 
in the late 1950's, which closely represent the median date associated 
wi'th the streamflow records in this study. 



Because the ability to accurately transfer streamflow data from 
gaged to ungaged sites is increased by raising network efficiencies, 
the application of remotely sensed land-use information to improve 
streamflow network models is a potentially valuable analytical tool. 
However, the generally favorable improvement in the network model of 
the Maryland district of the U.S. Geological Survey following inclusion 
of land-use data based on high-altitude photographs and satellite 
imagery may or may not be exceeded in other parts of the Nation. 
Accordingly, it is recommended that experiments, similar to those used 
in this report be conducted wherever remotely sensed land-use data are 
currently available. This would permit the making of accurate assess- 
ments of the use of remotely sensed land-use information to improve 
streamflow network models under a wide range of physiographic, 
climatic, and geologic settings. 
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Table A1 


T Streaioflow and basin characteristics. of stations used in multiple regression analysis, 
(see eKplanation on, p. 84), 
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, Table Al, — Screamflow and basin characteristics of stations used in multiple regression, analysis 
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249.9999 

289,9998 

12.9000 




Table A1 


Streaaflow and basin characteristics of stations used in aiultiple regression analysis 
— Continued. 


STATlOfJ 

/Jo 

Ca/ 

0 

0 

0 


0 


© 

0 

0 

~ 

0H64N6LT 

46HICJLT 

6A2F0«5T 



UDaoN.I 

uhuaniO 

FOkST.L 

FOrIST.8 

UftSN'.R^I 

l*778i10 

84.9000 

3.5000 

11.0000 

0.6000 

69.6000 

10.5000 

4.6000 

0.9003 

10. loo: 

50,3000 

l*7b0<!3 

20.9000 

59.9000 

19.2000 

0.0 

15.4000 

4.3000 

O.dOOO 

1.0003 

le.if oj 

3 D • 7 C Vv 

1476500 

3 k 0 V 00 

78.00011 

) 9.0000 

0.05G0 

1,7000 

i *0000 

0.3000 

C. 5000 

is.sdo: 

2,7000 

14e3?0(> 

c.v 

61 .dOOO 

37.O000 

0.3000 

O.U 

0.0 

0*0 

j * 0 

37.6000 

C.O 

UB350 0 

0.0 

52.3000 

17.7000 

0*0 

0.0 

0*0 

0.0 

0.2DDO 

17.500: 

0.0 

14943on 

0,9 

46.50nc 

52.5000 

0.6000 

0.0 

O.l) 

0.0 

o.roco 

51.8000 

0.0 

14B4500 

1.3000 

56.5000 

42. 2000 

0.0 

O.D 

1.3000 

O.D 

0.3 

42.2000 

1.3000 

1485000 

0.2C00 

49.4000 

50.2000 

0.0 

0.2000 

0.0 

0.0 

0,9000 

^9 . 3 0 0 V 

0.2000 

1465500 

0.2000 

19.2000 

79.6000 

0.0 

0.2000 

0.0 

0.0 

5.3000 

74.300C 

0,2000 

148600D 

o.c 

31.6000 


0.0 

0*0 

0.0 

o*c 

OtSODG 

67. »o:: 

0,0 

1466500 

5. 1000 

44, 3000 

49.6000 

0.400D 

l.vnoc 

3.2090 

0.0 

1.9009 

47 • 

5 » 1 G OO 

1467000 

1 .11)00 

57.4MC0 

41,3flJO 

u.o 

1 . OUCCf 

u.inoo 

c.o 

0.2003 

41.24c; 

1 .lOCC 

1487500 

0.0 

26.1000 

72. SO DO 

0.6U11Q 

0.0 

0.0 

0.0 

9 7005 

a5..cc: 

0.0 

146S500 

c.iooo 

88*0000 


0.0 

O.IOGO 

0.0 

0.0 

0 .2000 

4. 1 » ^ 0 V D 

0.1000 

1^69(^00 

OkO 

72.5UOO 

27.^njjn 

u.o 

0*0 

0*0 

0.0 

0*0 

27.5C0C 

0.0 

1490000 

0.0 

53.0000 

46*80DO 

0.2000 

0.0 

0.0 

G.O 

c.dOOO 

* Or j " 

0.0 

1491000 

0.1000 

55,80.)0 

43.9(it,o 

ft.o 

0,1000 

o.t) 

0.0 

0.9000 

A3.0D0D 

0,1000 

1492000 

0.0 

71.2000 

28.H000 

0.0 

0*0 

0*0 

0,0 

0.0 

28.6000 

0.0 

1492500 

0.0 

67.9000 

32.1000 

0.0 

0*0 

0.0 

G*0 

0.0 

32.1000 

c.o 

14V3000 

0.4000 

ro.iooo 

29.2000 

0.3000 

o.4(;o<; 

0.0 

0.0 

1.5000 

27.7000 

0,4000 

1493500 

0.0 

96.2030 

3tHr^oo 

0*0 

0*0 

0.0 

0.0 

0.0 

3.bC04^ 

5.0 

1494000 

0.0 

ra.bPiio 


a.p 

U.D 

0.0 

0.0 

o.aooo 

26.3..5C 

C*G 

1495C0O 

1.1000 

65,9000 

13.0000 

0.0 

0.9000 

0.2000 

O.D 

C.IOOO 

12.9C0V 

j • iDOU 

1495500 

1.2C00 

7V.b{><?0 

la.jfljo 

0,1000 

0. 700(1 

0.4000 

0.0 

3.1000 

15.2000 

I.IOOO 

1579000 

1.900U 

73.<»ono 

noo 

0,0 

1.9(100 

0*0 

Cl.O 

2.0COO 

22.700J 

1.9000 

1566000 

3.6000 

70.9000 

25,3000 

I). 1000 

2.3000 

1.2000 

0.1000 

1.6000 

23.7000 

3,5000 

1569390 

36.4000 

,;3. 7000 

40.0000 

0.0 

?7*30DO 

b. 7(*00 

2.4000 

6*7000 


33,0300 


D.D 

33.0000 

67,0000 

0.0 

0*0 

0.0 

0.0 

0*0 


0*0 

1591000 

1.50U0 

6&.30DO 

32.2000 


1.4000 

U.IOOO 

0*0 

(3.4000 

31.8DC0 

1.5000 

1 

la.ecoo 

38.6000 

42.91)1)0 

n.n 

12.6090 

3.8000 

2.1000 

1.5000 

41.4000 

16,4000 

1594600 

1.6000 

57.7000 

40.7000 

0.0 

1 . 600 (I 

0.0 

0.0 

0.0 

40,7500 

1.6000 

1 6&fi&0Q 

1.9000 

63.1000 

34.6000 

0*4000 

1 .6000 

0.3000 

0.() 

1*6000 

33.:300 

1.9C00 

1645200 

42.500(1 

40«VDO0 

16.6000 

0.0 . 

26.2001) 

14.00(10 

2.3000 

4.8000 

11.8000 

40.2000 

1658000 

7.2000 

24.7000 

66.0UU0 

0.1000 

3.50QD 


0.^000 

I.OOCO 

6T.COOC 

O.3000 

1653500 

63.2000 

4«80oa 


0*0 

39.7001' 

12.2000 

10.3000 

3.6000 

2d * Ao 0 5 

51.9000 

1 hft6550 

93.9000 

0.0 

6.1000 

0*0 

73.1)000 

12.0000 

9.0000 

3,7000 

2.4000 

, 55.0000 

1646000 

33.3CrOo 

25. 7000 

jooi; 

D.3000 

37.501)0 

10,0000 

5.8000 

2.6000 

iB.ioc; 

47.5000 

1649500 

45.1000 

15.8000 

38.9000 

0.2000 

24.90111) 

1 1.7000 

H.500D 

2.6030 

36.6CC9 

36.6000 

26S05U0 

Pb^DOOO 

42.4D0O 

31.6000 

0.0 

15,3000 

i.SOOO 

9.2000 

6.0000 

25.6000 

la. 8000 



Table A1 


, — Straamflow and basin characteristics 

— Continued. 



© 


© 

© 

Etttar 

A/O 

£RTSC^rK 

ERTStGHt 

ERTSF6ST 

EBTSw'iTS 

i?j47?eoo 

86.O000 

0.0 

14.0000 

C.O 

i47s:oo 

16.0000 

56,vOOC 

’ 30,3000 

2.3000 

IA7?E.'0 

I .0 

8C.0030 

20.30u0 

0.0 

lisrioc 

0 « 0 

51 .0000 

*9.0000 

C.O 

J 3S no 

c.c 

51.0000 

5.3000 

0.9 

i-SHSOD ’ 

0.7 

52,0000 

48.0000 

C.O 

j.si-soc 

I.mOOO 

n 1 . u 8 0 

3S • 0 0 0 0 

0.0 , 

1*55000 

0.0 

*3,0000 

56.0000 

1.0000 

1065500 

E . C 0 0 C 

2o . D 0 0 0 

72.0000 

8.0 

1066000 

0.0 

33.0030 

70.0000 

C.O 

1*66500 

4. COCO 

51,0000 

45.3000 

C 1 0 

14P7C00 

C.C 

5. .0000 

45.0000 

1 .9900 

14S7500 

0.0 

26. 0000 

72.0000 

0.0 

1*66500 

*.0000 

56,0000 

40.3000 

O.C 

j HcOOvO 

0.0 

71 .0000 

29.3000 

0.0 

14500GO 

C.G 

47.0800 

53.3000 

0.0 

1401000 

1,8000 

65.0000 

44.8000 

0.0 

I452Q00 

0.0 

53.0000 

7.9000 

C.O 

145ES00 

0.0 

65.0000 

31.9000 

0.0 

14S3000 

0.8 

74.0000 

26.0000 

0,0 

14S3500 

0.0 

97,0000 

3.3000 

0,0 

1*64000 

0.0 

71.0000 

29.9000 

0.0 

1455000 

C.O 

60.0000 

20.0000 

0.3 

1455500 

0.0 

62.0000 

36.9000 

0.0 

1575000 

0.0 

95.5000 

0,1000 

0.3 

1566000 

6.0000 

31 .0000 

1 7 . 0000 

0.3 

1585:00 

45,0000 

25 .0000 

30 .9000 

C.O 

15S0G00 

'5,0050 

13.0000 

<42.0000 

0.0 

1S51000 

0.0 

65.0000 

35,3000 

0.8 

155*500 

19.0C00 

35.0000 

46.0000 

O.C 

1554600 

l > . * 

58.0000 

42. 3000 

0.0 

J65S500 

0.0 

81, 9000 

19,0000 

0.0 

I64SZOO 

61. coco 

24.0000 

IS. 0000 

0.0 

16*6560 

54,0000 

0,0 

6.0000 

8.0 

164«»CaB 

49,0000 

34,0000 

16.0080 

1.0000 

1645SOO 

55.0000 

13.0000 

32.3000 

0.0 

1659500 

6.0000 

71.0000 

23.9000 

C.O 

16S3S00 

65,0000 

0.0 

15.0000 

0.0 

<7 1653C0C 

1.0000 

29.0000 

70,3000 

0.0 


stations used in multiple regression analysis 
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EXPLANATION 


B. 

Station No. 


Co]. 1 
Col. 2 

Col. 3 

Col. 4 
Col. 5 

Col. 6 
Col. 7 


Cols. 8-18 
Col. 19 


Col. 20 


Col. 21 


Col. 22 
Col. 23 


Missing data. 

These eight digit numbers are permanent nationwide 
numbers assigned by the U.S. Geological SurVfey fco 
stations at which streamflow data are collected on A 
recurrent basis. 

Drainage area, in square miles. 

Main-channel slope, in feet per mile, dete -mined from 
elevations at points 10 percent and 85 per ent of the 
distance along the channel from the gaging station to 
the drainage divide. 

Maxn-channel length, in miles, from the gaging station 
to the basin divide. 

Mean-basin elevation, in feet above mean sea level. 

Storage, in percent, of the drainage area covered by 
lakes, ponds, and swamps. 

Forest cover, in percent, of the drainage area covered 
by forests as shown on USGS 1:24,000 scale topographic 
maps. 

Soil index, a measure of potential maximum infiltration 
capacity, in inches, estimated from a map or from other 
data provided by the U.S. Soil Conservation_Serv.ice. 

Not used in the analysis. 

Mean annual precipitation, in inches, determined from 
an Isohyetal map prepared from National Weather Service 
records. 

Precipitation intensity, which is the maximum 24-hour 
rainfall, in inches, having a recurrence interval of 
2 years (24-hour 2-year rainfall) . 

Average annual snowfall, in inches, estimated from maps 
of average snowfall prepared from National Weather 
Service records. 

Average minimum January temperature, in degrees Fahrenheit. 
Average maximum July temperature, in degrees Fahrenheit. 
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Col. 24-28 

Col. 29 
Col. 30 
Col. 31-42 

i. 

Col. 43-54 
Col. 55-58 

f 

Col. 59-65 


Col. 66 

Col. 67 
Col. 68-71 

Col. 72-77 

Col. 78-81 


Flood-peak charac: eristics are represented by discharge 
from the annual f • ood-f requency curve at recurrence 
intervals of 2, 5, 10, 25, and 50 years. 

3 

Mean annual discharge, in ft /s. 

3 

Standard deviation of mean annual flows, in ft /s. 

3 

Mean monthly discharge, in ft /s beginning with Q-„ 
(October) . 

3 

Standard deviation on monthly flows, in ft /s. 

Low-flow characteristics are the annual minimum 7-day 
mean flows, in ft^/s at 2-year, 10-year, and 20-year 
recurrence intervals ( M iq» -q) > Col. 56 

not used. ’ * * 

Flood-volume characteristics represent the annual 
highest average flow, in ft^/s for 3-day periods at 
recurrence intervals of 2, 25, and 50 years and for 
7-day periods at recurrence intervals of 2, 10, 25, 
and 50 years. 

Fifty percentile discharge on the flow duration curve, 
in ft^/s. 

Not used in the analysis. 

Level I land use categories, in percent, determined 
from high altitude areal photographs. 

Level II land use categories, in percent, determined 
from high altitude areal photographs. , 

Level I land use categories, in percent, determined 
from Lands at (ERTS) imagery. 
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